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Abstract 

Probability  of  bit  ofijr  (P^)  performance  for  a  turbulent,  atmospheric,  optical 
communications  link  w4th  /convolutional  error  detection  and  correction  coding  was 
investigated.  The  cddes  investigated  were  all  rate  1/2  codes  with  constraint  lengths 
ranging  from  3  t  n^vas  assumed  that  the  coded  data  stream  was  sufficiently 
interleaved  s^that  the  turbulent  channel  could  be  considered  memoryless.  A  the¬ 
oretical  bound  on  Pi,  given  by  Viterbi  for  a  memoryless  channel,  was  used  to  plot 
curves  of  P/  versus  the  signal  intensity  of  the  laser  transmitter.  These  curves  were 
plotted  for  three  different  levels  of  atmospheric  turbulence  and  three  different  levels 


of  background  light. 
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USING  ERROR  DETECTION  AND  CORRECTION 
CODING  FOR  A  TURBULENT  ATMOSPHERIC 
OPTICAL  COMMUNICATIONS  LINK 


1.  Introduction 


1.1  Background 

John  Holt  made  it  clear  in  the  executive  S‘-*mmary  for  the  HAVE  LACE  test 
program  that  atmospheric  optical  commun'eations  hold  great  promise  for  the  Air 
Force. 


Laser  communications  offers  many  advantage.'  over  radio  frequency  (RF) 
systems  for  airborne  applications.  Most  of  these  advantages  stem  from 
the  frequency  of  the  carrier,  which  is  four  orders  of  magnitude  greater 
than  any  other  being  considered  for  airborne  use  The  higher  frequency 
allows  reduction  in  antenna  si'je  and/or  improvements  in  antenna  gain. 
High  gain  antennas  permit  these  systems  lo  pre  •  ide  long  range,  high  data 
rate  links  using  low  power  solid  state  transmitters.  Another  advantage  'A 
the  optical  frequency  is  the  ease  of  eliminating  the  backlobes  and  reducing 
the  sidelobes  by  orders  of  magnitude  o  ’er  RF  systems.  This  makes  laser 
communications  inherently  jam  resistant  and  difficult  to  intercept  (21). 


Large  data  rates,  low  probability  of  intercept  (LPl),  and  frequency  availability  are 
benefits  that  have  driven  the  fiber  optic.,  industry  through  the  roof,  but  only  recently 
hcis  a  lot  of  work  been  done  to  overcome,  or  least  work  around  the  problems  that 
atmospheric  effects  cause.  The  two  major  atmospheric  effects  are  the  extinction 
of  the  optical  beam  due  to  scattering  and  absorption,  and  the  disruption  of  the 
beam  due  to  clear-air  turbulence  (19).  The  Air  Force  test  program  HAVE  LACE 
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was  specifically  dimed  at  investigating  atmospheric  effects  on  high  altitude,  air- 
to-air  optical  communications.  Robert  Feldman  used  the  data  from  HAVE  LACE 
to  investigate  turbulence  induced  scintillation  of  the  optical  beam,  and  to  verify  a 
theoretical  model  for  describing  the  clear-air  turbulence  (13).  John  Holt  investigated 
the  effects  of  scattering  and  absorption  from  measurements  taken  in  the  HAVE  LACE 
test  flights  (20).  The  Air  Force  Geoph^  '‘cs  Lab  has  also  done  a  lot  of  work  in  the 
characterization  of  atmospheric  turbulence  (2,  4,  17,  46).  Hov/ever,  most  of  their 
experimental  work  to  date,  hjs  concentrated  on  the  ground- to-space  vertical  path. 
The  only  work  they  have  done  on  horizontal  paths  was  at  low  altitudes  (2,  34,  46). 
They  are  just  now  starting  to  study  the  use  of  a  laser  radar  for  remote  measurement 
of  turbulence  parameters  for  a  high  altitude,  horizontal  path  (3). 

In  order  to  predict  a  communications  systems’s  performance,  the  communica¬ 
tions  engineer  must  be  able  to  model  the  channel.  For  a  digital  laser  communications 
system,  predicting  performance  means  being  able  to  predict  the  probability  of  a  bit, 
or  symbol  error  at  the  receiver.  Modeling  the  channel  means  being  able  to  mathe¬ 
matically  characterize  atmospheric  turbulence  effects.  If  the  effects  of  the  channel 
on  the  signal  can  be  characterized,  various  techniques  are  available  to  combat  these 
effects. 

Feldman  suggested  in  his  thesis  that  error  detection  and  correction  (EDAC) 
coding  could  be  used  for  improving  performance  (13).  The  idea  for  this  thesis  was 
more  concisely  defined  through  discussions  with  Mr.  Feldman,  who  was  representing 
the  Wright  Research  and  Development  Center  (WRDC)  as  the  sponsor  for  this  thesis. 
A  typical  communications  system  is  designed  to  provide  a  certain  minimum  perfor¬ 
mance,  maximum  probability  of  error  in  the  digital  case,  under  the  worst  channel 
conditions  anticipated.  Improvement  in  the  channel  means  a  decreased  pro’oa’Dility 
of  error  at  the  same  range,  or  the  same  probability  of  error  at  an  increased  range.  In 
the  LPI  environment,  extended  range  or  improved  performance  is  not  necessarily  a 
benefit.  The  minimum  probability  of  error  required  for  reliable  decoding  of  the  signal 
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should  occur  at  the  receiver,  and  no  further.  Although  adaptive  coding  by  itself  does 
not  provide  a  LPI  capability,  when  combined  with  adaptive  power  control,  it  would 
allow  the  transmitter  to  be  operated  at  lower  power  levels  which  would  provide  a  LPI 
capability.  WRDC  wants  to  investigate  the  use  of  adaptive  ED  AC  to  combat  chang¬ 
ing  turbulence  in  the  atmospheric  channel.  If  the  turbulence  increases,  thus  raising 
the  probability  of  error,  then  the  EDAC  will  be  increased  to  provide  more  coding 
gain.  Conversely,  if  the  turbulence  decreases,  the  coding  gain  can  be  decreased  so  as 
to  prevent  the  signal  from  being  decoded  beyond  the  intended  receiver.  A  secondary 
benefit  of  adaptive  EDAC  is  the  maximization  of  data  throughput,  since  the  amount 
of  coding  redundancy  in  the  data  stream  is  minimized. 

1.2  Problem  Statement 

There  is  no  documented  study  of  the  use  of  adaptive  EDAC  in  an  atmospheric 
optical  communications  link  for  the  purpose  of  compensating  for  the  time  varying 
fading  caused  by  turbulence. 

1.3  Objective 

The  objective  of  this  thesis  is  to  show  that  EDAC  can  adaptively  compensate 
for  the  time  varying  fading  that  occurs  in  the  clear-air  turbulent  atmospheric  chan¬ 
nel,  allowing  the  communications  engineer  to  design  the  link  so  that  the  signal  is 
detectable  at  the  receiver,  but  not  beyond  the  receiver. 

1.4  Summary  of  Current  Knowledge 

Extensive  research  has  been  done  in  the  areas  of  atmospheric  turbulence  (1,2, 
4, 17, 23, 34,  41,  44),  and  optical  communications  in  space,  where  the  channel  is  con¬ 
sidered  memoryless.  In  a  memoryless  channel,  errors  occur  randomly  with  complete 
statistical  independence  between  them  (15,  36).  The  atmospheric  channel,  however, 
is  not  a  memoryless  channel.  Errors  occur  in  bursts  and  are  not  statistically  inde- 
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pendent.  This  makes  the  theoretical  analysis  of  the  probability  of  error  very  difficult 
(12:993).  As  a  result,  documented  research  into  the  use  of  EDAC  for  the  turbulent 
atmospheric  channel  is  sparse.  Most  of  it  is  experimental  in  nature,  comparing  ac¬ 
tual  measured  system  performance  to  the  theoretical  bounds  on  the  probability  of 
error.  These  bounds  are  derived  for  a  memoryless  channel  (5,  12,  25).  Cheng’s  dis¬ 
sertation  (6)  was  the  only  source  found  that  did  a  theoretical  performance  analysis 
of  an  optical  communications  system  over  a  turbulent,  atmospheric  channel.  He  did 
not  analyze  the  probability  of  error,  however.  His  measure  of  performance  was  the 
cutoff  rate  for  the  channel,  which  he  was  able  to  theoretically  determine. 

L5  Scope 

This  thesis  will  describe  a  communications  system  that  uses  a  direct  detection 
receiver.  Consequently,  the  turbulence  effect  of  interest  in  Chapter  II  is  scintillation 
since  it  affects  the  intensity  of  the  received  beam.  Phcise  distortions  will  not  need  to 
be  considered.  As  for  the  EDAC,  the  coding  gains  will  be  analyzed  for  convolutional 
codes  of  varying  constraint  lengths,  using  the  Viterbi  algorithm  (47:235-239)  for  soft 
decision  decoding.  The  constraint  length  of  the  convolutional  codes  will  be  limited 
to  10  since  present  technology  can  not  provide  real-time  decoding  with  the  Viterbi 
algorithim  for  anything  longer  (40:338). 

1.6  Methodology 

Due  to  the  “prohibitively  complex”  (12:993)  mathematics  needed  to  theoreti¬ 
cally  determine  the  prob.abilitj  of  bit  error  for  the  clear-air  turbulent  optical  channel, 
there  are  two  options  available  for  determining  the  probability  of  bit  error  curves. 
Experimental  measurements  can  be  taken,  or,  it  can  be  assumed  that  the  information 
signal  hcis  been  sufficiently  interleaved  so  as  to  make  the  channel  appear  memory¬ 
less.  Equipment  and  time  limitations  mandate  the  second  option  for  this  thesis. 
Davidson’s  results  (12)  will  be  used  to  determine  what  can  be  considered  sufficient 
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interleaving.  Davidson  showed  that  with  soft  decision  Viterbi  decoding  and  low  at¬ 
mospheric  turbulence,  an  interleaving  delay  of  about  5  milliseconds  was  sufficient  to 
bring  the  measured  probability  of  error  curves  veiy  close  to  the  curves  predicted  for 
a  memoryless  channel  (12:1001).  Assuming  a  memoryless  channel,  theoretical  values 
for  various  convolutional  codes  will  be  calculated  and  plotted.  Then,  assuming  suf¬ 
ficient  interleaving,  these  plots  can  be  used  to  compare  the  probability  of  bit  error 
for  the  different  coding  schemes. 

1.7  Materials  and  Equipment 

The  simplifications  resulting  from  the  assumptions  described  in  section  1.6  will 
allow  the  calculations  to  be  made  on  available  AFIT  computer  resources,  using  avail¬ 
able  software.  A  mathematical  analysis  package  such  as  MathCad,  Math  Station, 
or  Mathematica  will  be  needed,  as  well  as  a  matrix  oriented  analysis  tool  such  as 
MATRIXx  or  Matlab.  All  of  these  are  available  on  AFIT’s  Sun  workstations  except 
for  MathCad  which  is  available  on  the  AFIT  PC  network  server. 


1.8  Organization  of  the  Thesis 

Chapters  II  and  III  review  the  theory  needed  for  examining  a  digital  laser 
communications  system  that  uses  EDAC  to  combat  scintillation  effects.  Chapter  II 
covers  the  atmospheric  effects,  concentrating  on  turbulence  induced  scintillation, 
while  Chapter  III  covers  the  principles  of  EDAC,  concentrating  on  convolutional 
coding.  Chapter  IV  gives  a  detailed  description  of  the  methodology  used  to  derive 
the  equations  needed  for  calculating  the  theoretical  probability  of  bit  error  curves. 
Chapter  V  presents  the  results  of  the  calculations,  while  Chapter  VI  contains  the 
conclusions  and  recommendations  for  possible  follow  on  research  efforts. 
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11.  Turbulence  Theory 


2.1  Intr'^duction 

The  atmosphere  has  two  basic  effects  on  light  propagation,  extinction  and 
disruption  (19).  In  the  case  of  laser  communications,  extinction  means  that  energy 
from  the  laser  transmitter  is  prevented  from  reaching  the  receiver  or  detector,  while 
disruption  means  that  the  energy  that  reaches  the  detector  is  somehow  distorted  from 
the  form  it  had  when  it  left  the  transmitter.  Two  ways  that  the  atmosphere  can 
prevent  the  laser  energy  from  reaching  the  detector  are  absorption  and  scattering. 
Water,  carbon  dioxide  and  other  molecules  in  the  atmosphere  will  absorb  light, 
although  the  amount  of  absorption  is  wavelength  dependent.  Figure  2.1  shows  that 
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Figure  2.1.  Absorbtion  Spectrum 

there  are  certain  wavelength  regions  where  absorption  is  minimized.  These  regions 
are  known  as  absorption  windows.  The  proper  selection  of  wavelength  will  allow 
the  system  to  operate  in  one  of  these  windows.  Particles  in  the  atmosphere,  called 
aerosols,  will  cause  scattering  of  the  laser  energy,  preventing  any  mecisureable  amount 
of  energy  from  reaching  the  detector.  If  absorption  and  scattering  can  be  considered 
negligible,  then  the  atmosphere  is  described  as  being  “clear  air  turbulent”  (CAT) 
(18).  Goodman  does  a  good  job  of  distinguishing  between  the  CAT  atmosphere,  and 
an  atmosphere  where  absorbtion  and  scattering  are  not  negligible. 
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We  shall  always  assume  that  the  scale  sizes  are  much  larger  than  the 
wavelength  of  the  radiation  being  used.  This  assumption  eliminates  from 
consideration  problems  involving  imaging  through  clouds  or  aerosols,  for 
which  the  scale  sizes  of  the  inhomogeneities  are  comparable  with  or 
smaller  than  an  optical  wavelength  and  for  which  the  refractive  index 
changes  are  sharp  and  abrupt.  This  latter  class  of  problems  may  be  re¬ 
ferred  to  as  “imaging  through  turbid  media”  whereas  we  are  concerned 
here  with  “imaging  through  turbulent  media”  for  which  the  refractive 
index  changes  are  smoother  and  coarser  (18:362). 


The  “sharp  and  abrupt”  changes  in  the  index  of  refraction  cause  large  changes  in 
the  direction  of  the  radiation,  which  is  nothing  more  than  scattering.  Goodman 
also  makes  the  assumption  that  absorption  is  negligible  in  his  analysis  of  the  CAT 
atmosphere  (18:385). 


S.2  Turbulence  Effects  on  Laser  Beam  Propagation 

Knowing  what  conditions  are  necessary  to  define  a  CAT  atmosphere,  what  is 
the  effect  of  such  an  atmosphere  on  the  propagation  of  the  laser  beam?  The  two  pri¬ 
mary  effects  are  scintillation  and  coherence  degradation.  Pratt  defines  scintillation 
as  “small  scale  destructive  interference  within  the  beam  cross-section,  causing  vari¬ 
ations  in  the  spatial  power  density  at  the  receiver  (37:133).”  This  same  interference 
also  causes  a  coherence  degradation  since  the  interference  is  happening  randomly. 
However,  in  this  thesis,  the  use  of  a  direct  detection  receiver  is  assumed  which  is 
usually  the  case  in  present  laser  communications  links.  Since  only  intensity  varia¬ 
tions  affect  the  probability  of  bit  error  with  this  type  of  receiver,  coherence,  or  phase 
degradation  of  the  beam  is  not  a  problem. 


2.8  Modeling  of  Scintillation 

Propagation  of  light  in  the  atmosphere,  just  as  everywhere  else,  can  be  de¬ 
scribed  by  Maxwell’s  equations.  Assuming  that  the  atmosphere  heis  zero  conductiv¬ 
ity  and  unit  magnetic  permeability,  and  that  the  electric  field  is  sinusoidal  with  a 
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time  dependence  of  Maxwell’s  equations  reduce  to 


V-H 

=  0 

(2.1) 

V  X  E 

=  ikH 

(2.2) 

WxH 

=  -ikn^E 

(2.3) 

V  •  {n'^E) 

=  0 

(2.4) 

where  E  and  H  are  the  electric  and  magnetic  fields  respectively,  n  is  the  index  of 
refraction,  k  is  the  wave  number,  w/c,  where  a;  is  the  radian  frequency  and  c  is  the 
speed  of  light,  and  i  is  y/^.  Using  vector  calculus,  these  equations  can  be  reduced 
to  the  more  familiar  wave  equation: 


V^E  +  k‘^n‘^E  +  2V{E-Vlogn)=0  (2.5) 

The  Icist  term  on  the  left  hand  side  describes  the  depolarization  effects.  This  term 
can  be  considered  negligible  for  wavelengths  (A)  <  lo  where  lo  is  the  inner  scale  of 
the  turbulence  (10,  27,  38,  42,  43,  45).  For  atmospheric  turbulence,  /„  >  1mm,  so 
with  optical  waves  where  A  is  on  the  order  of  a  micron,  this  is  a  valid  assumption. 
If  you  also  consider  the  components  of  the  electric  field  separately,  equation  2.5  can 
be  written  in  its  scalar  form  as: 


V'^E  +  k^n^E  =  0  (2.6) 

If  the  channel  was  the  space  channel,  where  n  «  1,  or  in  a  fiber  with  n  «  1.55, 
then  solving  this  equation  is  fairly  straightforward.  Unfortunately,  the  atmosphere, 
due  to  turbulence,  is  inhomogeneous  and  anisotropic.  As  a  result  n  is  neither  con¬ 
stant,  nor  even  deterministic.  As  Clifford  states: 


The  refractive  index  fluctuations  that  we  are  considering  are  those  that 
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result  from  the  naturally  occurring  random  fluctuations  in  wind  velocity 
called  turbulence.  Because  turbulence  is  a  random  process,  it  must  be 
described  in  terms  of  statistical  properties.  It  is  hopeless  to  attempt  to 
predict  the  velocity  of  a  given  parcel  of  air  as  it  is  buffeted  about  in 
the  open  atmosphere.  In  fact,  it  is  nearly  impossible  to  do  a  complete 
statistical  analysis.  Therefore,  we  settle  for  a  more  meager  description. 

(9:11) 

Obviously  the  first  tcisk  is  to  statistically  characterize  n.  The  usual  way  of 
describing  a  random  process  is  in  terms  of  its  mean  and  autocorrelation  functions 
(44:3).  Assuming  wide  sense  stationarity,  the  Fourier  transform  of  the  autocorrela¬ 
tion  function  can  then  be  used  to  describe  the  power  spectral  density  of  the  process. 
A  good  example  of  this  is  the  characterization  of  white  noise.  If  white  noise  is  sam¬ 
pled  at  any  two  time  instants,  there  is  complete  statistical  independence  between 
the  samples.  Therefore,  the  autocorrelation  is  zero.  Only  when  the  time  difference 
between  samples,  t,  is  zero,  is  the  autocorrelation  non-zero.  The  autocorrelation 
function  is  then  given  by  a  delta  function  at  r  =  0  with  a  magnitude  that  is  related 
to  the  strength  of  the  fluctuations  in  the  noise,  given  by  the  variance  of  the  gaussian 
noise,  cis  shown  in  Figure  2.2.  The  Fourier  transform  of  this  delta  function  is  a  flat, 

Rxxi,'^) 


- 1 -  r 

0 

Figure  2.2.  Autocorrelation  Function  for  White  Noise 


infinitely  wide  power  spectral  density  whose  magnitude  is  inversely  proportional  to 
the  strength  of  the  fluctuations  in  the  noise.  Th*^  power  spectral  density  is  shown 
in  Figure  2.3.  As  was  mentioned  above,  one  condition  for  being  able  to  take  the 
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Figure  2.3.  Noise  Power  Spectral  Density 

Fourier  transform  is  that  the  random  process  be  wide  sense  stationary.  That  is,  its 
mean  value  must  be  independent  of  time,  and  its  autocorrelation  must  be  a  function 
only  of  the  time  difference,  t  =  t\~  t2,  not  the  actual  times  U  and  ^2- 

Just  as  the  power  spectral  density  of  white  noise  is  used  in  analyzing  RF 
communications  links,  a  power  spectral  density  for  the  fluctuations  in  n  is  needed 
in  order  to  analyze  the  atmospheric  optical  communications  link.  Unfortunately, 
most  real  random  processes  are  not  stationary,  so  as  in  all  applications  of  theory 
to  the  real  world,  the  question  is,  how  good  is  the  approximation?  Herein  lies  the 
problem  in  characterizing  the  random  fluctuations  in  n.  It  is  not  even  close  to  being 
a  stationary  random  process. 

Atmospheric  turbulence  always  contains  large  scale  components  which 
usually  destroy  the  homogeneity  and  isotropy  of  fields  of  the  meterological 
variables;  moreover,  these  components  cause  the  meteorological  fields  to 
be  non-stationary.  (44:19) 

In  general,  the  random  fluctuations  of  n  are  a  function  of  position,  wavelength, 
and  time:  n  =  n(f,  A,f).  However,  the  wavelength  dependence  is  usually  negligible 
in  the  atmosphere.  The  time  dependence  is  also  negligible  since  “The  time  required 
for  light  to  propagate  through  the  atmosphere  is  only  a  fraction  of  the  fluctuation 
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time  of  the  random  refractive  index  component.”  (18:385-6)  All  that  is  left  is  the 
spatial  variation  of  n,  n(f).  The  autocorrelation  function  of  n  can  be  written  as: 


5„(r)=£;{n(r1)n(r2)}  (2.7) 

where  f=fl  — r2.  If  n  was  a  stationary  process,  then  the  power  spectral  density  of  the 
fluctuations  could  be  obtained  by  Fourier  transforming  equation  2.7.  Unfortunately, 
as  mentioned  earlier,  n  is  not  a  stationary  process.  The  mean  value  changes  with 
r,  and,  due  to  the  inhomogeneities  of  the  atmosphere,  the  autocorrelation  function 
depends  not  only  on  the  vector  distance  between  sample  points,  rl  —f^,  but  on  the 
location  of  the  sample  points,  r1  and  r2,  as  well. 

Kolmogorov’s  solution  to  this  problem  was  the  structure  function.  The  struc¬ 
ture  function  for  n  is  the  square  of  the  average  difference  in  n  at  the  two  points, 
rather  than  the  average  of  the  product  of  n  at  the  two  points  as  is  the  autocorrela¬ 
tion  function.  Tatarski  describes  the  benefit  of  the  structure  function  cis: 


The  difference  between  the  values  of  the  /(f)  at  two  points  fj  and  r2 
is  chiefly  affecte-^  only  by  inhomogeneities  of  the  field  f  with  dimensions 
which  do  not  exceed  the  distance  |fl  —  f^j.  If  this  distance  is  not  too 
large,  the  largest  inhomogeneities  have  no  effect  on  /(fi)  —  /(fa)  and 
therefore  the  structure  function 

I>/(n,r5)=[/(r;)-/(r;)]" 

can  depend  only  on  fj  —  f2.  At  the  same  time  the  value  of  the  autocorre¬ 
lation  function  is  affected  by  inhomogeneities  of  all  scales,  so  that  for  the 
same  values  of  fi  and  f2  the  correlation  function  can  depend  on  each  of 
the  arguments  separately  and  not  just  on  the  difference  fj  —  f^.  (44:19) 

The  structure  function  takes  advantage  of  a  process  which  has  stationary  increments. 
This  is  true  for  n  since  the  atmosphere  can  be  considered  locally  homogeneous  within 
the  distance  |fj  —  f^j.  In  describing  turbulence,  whether  it  is  atmospheric,  oceanic  or 
some  other  medium,  the  parameter  whose  fluctuations  are  most  apparent  is  velocity. 
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In  an  airplane,  the  velocity  changes  are  felt  as  buffeting.  In  the  water,  currents 
and  eddies  display  the  fluctuating  nature  of  the  velocity.  Other  parameters  that 
depend  directly  or  indirectly  on  the  velocity  changes,  can  be  understood  better  by 
understanding  the  velocity  fluctuations.  Kolmogorov  showed  that  as  long  eis  —  7^| 
is  small  enough,  the  structure  function  of  the  velocity  is  given  by 

Drr  =  Cyf^  (2.8) 

where  is  defined  as  the  structure  constant  of  the  velocity.  The  structure  constant 
is  a  measure  of  the  strength  of  the  fluctuations  (9:12),  just  as  the  variance  of  gaussian 
noise  is  a  measure  of  the  strength  of  the  fluctuations  in  the  noise.  Equation  2.8  is 
similar  to  Figure  2.2  in  that  Drr  is  proportional  to  C^,  while  Rxx{t)  is  proportional 
to 

Unfortunately  there  is  no  direct  relationship  between  the  velocity  fluctuations 
and  the  fluctuations  in  n.  The  index  of  refraction  depends  on  pressure,  temperature, 
and  wavelength  as  shown  in  equation  2.9  (9:10). 

n  -  1  =  77.6(1  +  7.52  •  10-^X-^){P/T)  ■  10"°  (2.9) 

where  P  is  given  in  millibars,  T  is  in  degree.^  Kelvin  and  A  is  in  microns.  The  left 
hand  side  of  the  equation  is  written  in  a  form  which  directly  gives  the  value  of  the 
fluctuation  of  n  eis  it  varies  about  its  mean  free  space  value  of  1.  Equation  2.9  is  only 
good  for  over  land  paths  since  it  neglects  the  effects  of  humidity.  Over  the  water,  a 
pressure  humidity  product  factor  makes  an  appreciable  contribution  (9). 

Corrsin  showed  that  potential  temperature  obeys  the  same  two-thirds  law  as 
the  velocity  fluctuations  (11).  How  then  can  n  be  related  to  potential  temperature? 
Clifford  describes  the  quantity  potential  temperature  as:  “. . .  the  difference  between 
the  absolute  temperature  T  and  the  changes  in  temperature  with  height  caused  by 
an  adiabatic  lapse  rate  ...0  =  T  —  7o2”  (9:12).  Taking  the  partial  derivative  of 
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equation  2.9  with  respect  to  pressure  and  temperature  and  evaluating  at  A  =  .6/xm, 
which  corresponds  to  red  light  (9:13)  we  get: 

6n  =  19P/T  {SP/P  -  8T/T)  •  10"®  (2.10) 

If  it  is  also  assumed  that  the  pressure  variations  are  negligible  for  a  fixed  altitude 
then  the  fluctuations  in  n,  at  a  given  height,  are  a  function  only  of  the  fluctuations  in 
temperature  at  that  height,  and  the  assumed  constant  pressure  at  that  height.  At  a 
fixed  altitude,  the  fluctuations  in  actual  temperature  also  obey  the  |  law.  Therefore 
6n  does  also  (9:13).  Consequently,  the  structure  function  for  n  can  be  written  as: 

D„(r)  =  Cyl^  (2.11) 

where  is  given  by  (9:13): 

01  =  X  (2.12) 

Although  the  fluctuations  of  n  itself  are  not  a  stationary  random  process,  the  fluc¬ 
tuations  of  n  about  its  mean  value,  which  from  here  on  will  be  designated  cis  n',  is  a 
stationary  process  if  the  atmosphere  can  be  considered  locally  homogeneous  within 
the  same  constraints  on  —  r^l  that  applied  to  the  structure  function.  The  power 
spectral  density  of  the  fluctuations  can  now  be  found  by  Fourier  transforming  the 
autocorrelation  function  of  n'  (9:18). 

*»'(«)  =  (2^  /  (2.13) 

Assuming  that  the  atmosphere  can  be  considered  locally  isotropic.  Equation  2.13 
can  be  rewritten  in  scalar  form  as  : 

1  t°° 

^n'(«^)=2^yo  ‘>'Bn'{r)sin{Kr)dr  (2.14) 
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Again  using  the  assumption  that  the  atmosphere  is  locally  homogeneous  and  isotropic, 
equation  2.14  can  be  written  in  terms  of  the  structure  function  Dn'{r)i  since  under 
these  conditions: 

B„»(0)  -  Bn'{r)  =  ~Dn'{r)  (2.15) 

Combining  equations  2.14,  and  2.15,  Dn'{r)  can  be  written  in  terms  of  as: 


TOO 

Dn'{r)  =  StT  / 

Jo 


1- 


sm(/cr) 


Kr 


dK 


(2.16) 


Inverting  equation  2.16  is  not  eis  trivial  as  Fourier  transforming  a  delta  function,  but 
it  can  be  done  with  the  judicious  use  of  derivatives  (44)  giving: 


_  ,  .  1  sin(Kr)  d  , 


r^^^DAr) 


dr 


(2.17) 


Substituting  equation  2.11  into  equation  2.17,  and  integrating,  $n'(K)  can  be  written 
in  terms  of  C^;: 

$„/(«)  =  (2.18) 

Equation  2.18  is  governed  by  the  conditions  of  local  homogeneity  and  isotropy.  The 
resulting  limits  on  k  are  given  by  Kolmogorov’s  inertial  subrange  model,  27r/I,o  ^ 
K  <C  27r//o,  where  lo  is  the  inner  scale  of  the  turbulence  and  Lo  is  the  outer  scale 
of  the  turbulence.  Clifford  states  that  :  “Near  the  ground  lo  is  on  the  order  of 
millimeters  and  Lo  of  the  order  of  meters.  Both  quantities  appear  to  increase  with 
height  above  the  ground  (9:12).”  Equation  2.18  needed  some  corrections  in  order  to 
accurately  model  the  physical  reality  of  the  power  spectral  density  of  the  fluctuations 
in  n'.  Goodman  descibes  these  corrections  cis:  (18:389) 


When  /c  reaches  ...  the  form  of  again  changes.  Turbulent 

eddies  smaller  than  a  certain  scale  size  dissipate  their  energy  as  a  result 
of  viscous  forces,  resulting  in  a  rapid  drop  in  for  k  >  /c„,.  The 
scale  size  —  27r//Cm  is  referred  to  as  the  inner  scale  of  the  turbulence. 
A  typical  value  for  L  near  the  ground  is  a  few  millimeters.  Tatarski 
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includes  the  rapid  decay  of  #n'('')  >  «m  by  use  of  the  model 

#n'(«)  =  (2.19) 


This  equation  is  a  reasonable  approximation  provided  Km  is  chosen  to 
equal  5.92//o  and  k  >  «<,• 

Spectra  2.18  and  2.19  both  have  nonintegrable  poles  at  the  origin.  In 
fact,  since  there  is  a  finite  amount  of  air  associated  with  the  Earth’s 
atmosphere,  the  spectrum  can  not  become  arbitrarily  large  as  k  0. 
To  overcome  this  defect  in  the  model,  a  form  known  as  the  von  Karman 
spectrum  is  often  adopted.  The  spectrum  is  then  expressed  as 


$„»(«:) 


0.033C2,  4- 

(/c2  + 


(2.20) 


A  plot  of  equation  2.20,  showing  the  inertial  subrange  limits  is  given  iu  Figure  2.4. 
4>„»,  normalized  by  C^.,  is  plotted  against  the  spatial  frequency,  k,  which  is  inversely 
proportional  to  a  spatial  distance. 


K 


Figure  2.4.  Power  Spectral  Density  of  Fluctuations  in  n 

(18:390) 
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Clifford  then  goes  through  the  solution  of  equation  2.6  using  the  .  ?thod  of 
small  pertubations  and  Rytov’s  approach  of  taking  the  log  of  the  ampl,  ie  and 
representing  it  as  x  (9:26-34).  The  resulting  equations  for  the  power  spectral  densities 
of  the  log-normal  amplitude,  ai  1  phase  spectra,  F,  are  given  in  equation  2.21. 


0) 

Fs(/c,  0) 


=  irk^L 


(2.21) 


The  minus  sign  gives  the  two  dimensional  power  spectral  density  of  the  log-normal 
amplitude  spectra  and  the  plus  sign  gives  the  two  dimensional  phase  spectra.  The 
justifications  for  Rytov’s  transformation  is  simply  that  the  resulting  log-normal  dis¬ 
tribution  of  the  amplitudes  matches  the  empirical  data  for  low  levels  of  turbulence. 
As  the  level  of  turbulence  increases,  either  because  of  increasing  C^i  or  longer  path 
length,  this  model  breaks  down  as  the  received  signal  can  no  longer  be  described 
by  a  log-normal  distribution.  More  research  has  been  done  to  extend  the  model 
(42:45-106),  but  th't  is  beyond  the  scope  of  this  thesis. 
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III.  Error  Detection  and  Correction  Coding 


3.1  Introduction 

EDAC  got  its  real  start  with  Shannon’s  theorem  on  channel  capacity  in  1948. 
Shannon  showed  that  for  a  data  rate  that  was  less  than  the  capacity  of  the  channel, 
an  arbitrarily  small  probability  of  bit  error  could  be  acheived  through  the  use  of 
EDAC. 

Let  a  discrete  channel  have  the  capacity  C  and  a  discrete  source  the 
entropy  per  second  H.  If  H  <  C  there  exists  a  coding  system  such  that 
the  output  of  the  source  can  be  transmitted  over  the  channel  with  an 
arbitrarily  small  frequency  of  errors  (39:39). 

Shannon’s  theorem  said  nothing  about  how  to  do  the  coding,  just  that  it  could  be 
done.  The  development  of  the  codes  Wcis  left  to  others.  There  are  several  good 
textbooks  on  the  principles  of  EDAC  (8,  28,  40,  47).  This  chapter  will  provide  a 
summary  of  the  theory  needed  for  understanding  the  use  of  EDAC  in  a  clear-air 
turbulent,  atmospheric  optical  communications  link. 

3.2  Block  Codes  vs.  Convolutional  Codes 

There  are  two  basic  types  of  codes  that  are  used  in  existing  communications 
systems,  block  codes  and  convolutional  codes.  According  to  Viterbi,  convolutional 
codes  “invariably  outperform  block  codes  of  the  same  order  of  complexity”  (49). 
Convolutional  codes  are  also  easier  to  decode  with  a  soft  decision  decoder  than 
block  codes  (40  ’31).  Soft  decision  decoding  is  simply  a  means  of  obtaining  more 
information  from  the  output  of  the  demodulator.  Hard  decision  decoding  puts  out  a 
one  or  a  zero  based  on  the  relationship  between  the  output  level  and  the  threshhold 
level.  If  the  output  is  above  the  threshhold,  a  one  is  output.  If  the  output  is  below 
the  threshhold  a  zero  is  output.  No  information  is  obtained  about  how  “good”  the 
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decision  is.  Soft  decision  decoding  quantizes  the  output.  Instead  of  only  having  one 
threshhold,  there  are  2^'  —  1  threshholds,  where  k  is  the  number  of  quantization  bits. 
The  output  can  now  contain  information  on  how  close  the  output  is  to  a  zero  or 
a  one.  Figure  3.1  demonstrates  this  for  ^  =  3  (40:331).  Eight  level  quantization 


000  001  010  on  100  loi  no  in 

Figure  3.1.  Soft  Decision  Decoding 

{k  =  3)  is  usually  used  since  it  gives  a  2  dB  improvement  in  performance,  whereas 
infinite  quantization  only  provides  another  .2  dB  of  improvement  (40:331).  Soft 
decision  decoding  has  been  shown  to  provide  between  1  and  2  dB  more  coding  gain 
over  hard  decision  decoding,  when  used  on  a  turbulent  atmospheric  channel  (12). 
These  factors  have  driven  the  use  of  convolutional  codes,  with  soft  decision  decoding, 
for  use  on  an  atmospheric  optical  communications  link. 

3.3  Convolutional  Code  Description  Methods 

The  usual  implementation  of  a  convolutional  encoder  is  a  sequence  of  shift 
registers  with  modulo-2  adders  connected  to  certain  registers  in  the  sequence.  A 
classic  example  used  in  several  textbooks  (8,  40,  47),  was  first  used  by  Viterbi  (49). 
It  is  shown  in  Figure  3.2.  The  code,  or  encoder,  is  described  by  representing  the 
taps  octally.  This  code  for  example,  would  be  represented  by  (5,7),  since  the  top 
modulo-2  adder  has  taps  coming  from  the  third  bit,  (numbered  from  right  to  left) 
and  from  the  first  bit.  In  binary  this  would  be  101,  octally  a  5. 
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Figure  3.2.  K  =  3  Encoder  Example 


The  code  can  also  uo  represented  by  a  tree,  as  shown  in  Figure  3.3.  Usually 
the  encoder  is  assumed  to  start  with  all  the  shift  registers  containing  zeroes.  An 
input  bit  of  zero  is  represented  by  the  upper  branch  while  a  one  follows  the  lower 
branch.  The  output  codewords  are  shown  on  each  branch.  As  an  example,  an  input 
message  of  10101  would  result  .n  a  transmitted  message  of  0001101011  where  the 
rightmost  bit  of  the  input  message  is  the  first  bit  into  the  encoder  and  the  rightmost 
bit  of  the  transmitted  message  is  the  first  bit  out.  The  path  for  this  message  is 
shown  in  Figure  3.3.  The  tree  diagram  becomes  unwieldy  for  an  input  sequence  of 
any  appreciable  length,  since  the  number  of  branches  increases  as  2^,  where  L  is  the 
length  of  the  input  message.  Fortunately,  however,  the  tree  starts  to  repeat  itself 
after  K  branchings,  where  K  is  the  number  of  shift  registers  in  the  encoder,  known  as 
the  constraint  length  of  the  code.  As  can  be  seen  in  Figure  3.3,  the  fourth  branching 
duplicates  the  eight  codewords  from  the  upper  and  lower  halves  of  the  tree.  This 
means  that  the  output  codeword  doesn’t  depend  on  the  K+1  previous  bit.  Thus  the 
constraint  length  determines  the  memory  of  the  encoder.  This  feature  allows  the 
tree  diagram  to  be  simplified  to  a  trellis  diagram  as  shown  in  Figure  3.4.  The  trellis 
diagram  shows  the  possible  state  transitions,  the  resulting  codeword,  and  the  input 
message  bit  that  caused  that  state  transition.  An  input  bit  of  0  is  represented  by 
a  solid  line,  and  an  input  bit  of  1  is  reprented  by  a  dcished  line.  As  an  example, 
again  assuming  the  encoder  is  started  with  all  zeroes,  the  two  possible  transitions 


Figure  3.3.  Tree  Diagr.im  for  Example  Encoder 
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Figure  3.4.  Trellis  for  Example  Encoder 

are  to  state  10,  or  to  remain  in  state  00.  The  transition  to  10,  caused  by  an  input 
bit  of  one,  results  in  an  output  codeword  of  11,  whereas  an  input  bit  of  0  causes  the 
encoder  to  remain  at  the  00  state,  resulting  in  an  output  codeword  of  00.  It  will 
be  discussed  later  how  useful  the  trellis  diagram  is  in  visualizing  the  decoding  of  a 
convolutional  code,  but  for  now  a  clearer  picture  of  the  state  transitions  is  given  by  a 
state  diagram.  The  state  diagram  for  the  encoder  of  Figure  3.2  is  given  in  Figure  3.5. 

The  state  diagram  is  the  most  compact  way  of  representing  the  encoder  since 
its  size  depends  only  on  the  number  of  states  and  not  on  the  length  of  the  input 
message  as  the  tree  and  trellis  diagrams  do.  The  number  of  states  is  given  by  2^“^ 
so  it  is  still  going  to  be  complicated  for  K  >  i.  A  solid  line  still  represents  an 
input  bit  of  0,  and  the  dashed  line  still  represents  a  1,  as  in  the  trellis  diagram.  All 
possible  state  transitions  are  shown  with  their  respective  output  codewords.  The 
diagrams  shown  in  figures  3.3,  3.4,  and  3.5  are  patterned  after  those  given  by  Sklar 
(40:322-326).  The  only  difference  is,  the  taps  on  Sklar’s  encoder  are  given  octally 
by  (7,5)  instead  of  (5,7),  causing  the  order  of  the  bits  in  the  output  codeword  to  be 
swapped.  The  codewords  00  and  11  are  thus  unchanged,  while  the  01  becomes  a  10, 
and  vice  versa. 
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00 


01 

Figure  3.5.  State  Diagram  for  Example  Encoder 


As  mentioned  earlier,  the  trellis  is  a  useful  diagram  for  visualizing  the  decod¬ 
ing  process.  This  is  demonstrated  in  Figure  3.6  which  shows  both  the  encoder  and 
decoder.  The  branches  of  the  encoder  trellis  diagram  are  still  labeled  with  the  out¬ 
put  codeword  resulting  from  that  state  transition  path.  The  decoder  trellis  branches, 
however  are  labeled  with  an  appropriate  path  metric  describing  the  conditional  prob¬ 
ability  of  that  path  being  the  correct  path  given  the  input  codeword.  The  bcisis  of 
Viterbi  decoding  can  be  summed  up  in  one  sentence.  “If  any  two  paths  in  the  trellis 
merge  to  a  single  state,  one  of  them  can  always  be  eliminated  in  the  search  for  the 
optimum  path  (40:334).”  The  elimination  is  done  by  comparing  the  cumulative  path 
metric  for  the  two  paths.  The  path  with  the  higher  metric  will  be  discarded.  If  the 
path  metrics  are  equal,  one  of  them  is  discarded  by  an  arbitrary  rule.  The  Viterbi 
algorithm  has  been  shown  to  be  an  optimum  decoding  scheme  (35).  This  procedure 
can  be  seen  at  time  <4  of  the  decoder  trellis  in  Figure  3.6.  In  this  example  the  Ham¬ 
ming  distance,  defined  as  the  difference  in  the  number  of  ones  between  the  received 
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Encoder  Trellis 


Input  data :  j  1  0  1  1 

Output:  11  10  10  00  10 

Received :  11  10  10  01  10 

a  =  00  1^  2 

0\ 

b  =  10  + 

c  =  01  ^ 

d  =  ll  + 

Decoder  Trellis 


Figure  3.6.  Encoder  and  Decoder  Trellis  Diagram 
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codeword  and  the  transmitted  codeword,  is  used  as  the  path  metric.  Typically  in  a 
optimal  decoder  the  metric  that  is  used  is  some  type  of  maximum  likelihood  metric, 
but  for  simplicity’s  sake,  this  example  uses  the  Hamming  metric.  There  are  two 
paths  merging  at  state  a,  at  time  U.  The  first  path,  which  stayed  in  state  a  for 
three  time  periods,  has  a  cumulative  path  metric  of  4.  The  other  path  which  goes 
through  states  b  and  c  before  returning  to  state  a,  heis  a  cumulative  path  metric  of 
3.  The  first  path  has  a  higher  path  metric,  and  therefore  can  be  eliminated.  As  the 
decoder  progresses  deeper  into  the  trellis,  more  paths  are  eliminated  until  only  one 
path  remains  from  time  ti  to  time  t2,  and  the  first  bit  can  then  be  output.  Figure  3.7 
shows  the  surviving  paths  after  three  and  four  codewords  have  been  received.  The 
cumulative  path  metrics  for  the  upper  and  lower  paths  arriving  at  a  node  are  shown 
to  the  right  of  the  final  node.  For  example,  at  time  tg  the  node  for  state  a  has  an 
upper  path  whose  cumulative  path  metric  is  4  and  a  lower  path  whose  cumulative 
path  metric  is  1.  Figure  3.8  shows  the  surviving  paths  at  times  tg  and  h.  Since  there 
are  still  two  paths  going  from  time  ^2  to  time  is,  no  output  bit  can  be  decided  yet. 
Figure  3.9  shows  the  surviving  paths  at  time  ts.  The  codeword  received  at  time  ts 
allows  the  elimination  of  all  but  one  path  between  time  and  <4,  allowing  the  output 
of  a  zero  for  the  third  message  bit.  Figure  3.10  shows  that  no  more  output  bits  can 
be  decided  upon,  since  there  are  still  two  paths  between  time  ^4  and  ts.  Figure  3.11 
demonstrates  that  there  still  is  not  enough  progression  into  the  trellis  to  decide  on 
another  output  bit.  Figure  3.12  shows  that  more  progression  is  still  needed.  Fig¬ 
ure  3.13  shows  that  after  picking  the  paths  with  the  higher  cumulative  path  metrics 
at  time  <11,  only  one  path  remains  for  the  five  bit  durations  that  correspond  to  the 
five  message  bits.  It  also  shows  that  the  error  in  the  fourth  codeword  was  corrected. 
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message  bit:  110 

codeword  sent:  11  10  10 


output  bit:  "I  ‘I  "I 

message  bit:  110  1 

codeword  sent:  11  10  10  00 

codeword  received:  11  10  10  01 


Figure  3.7.  Surviving  Paths  at  <4,^5 
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message  bit: 
codeword  sent: 
codeword  received: 

+ 


+ 


+ 

output  bit: 


+ 


+ 


+ 

output  bit: 


Figure  3.8.  Surviving  Paths  at  <6,^7 
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COlKd  coico  coico 


Figure  3.9.  Surviving  Paths  at  is 


message  bit:  110110  00 

codeword  tx:  11  10  10  00  10  10  01  01 

codeword  rx:  11  10  10  01  10  10  01  01 

K 

\ 

+ 

+ 

+ 

output  bit:  110???' 

Figure  3,10.  Surviving  Paths  at  <9 
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message  bit:  110  110  0 

codeword  sent:  11  10  10  00  10  10  01 

codeword  received:  11  10  10  01  10  10  01 

%  h 

\  ^  / 

+  X  +  +  /a  +  X/  ^ 

+  +  \^  +  /  +  \^  +  /  N 

+  +  '"-K  +  +  'X  +  4 

output  bit:  11????? 


Figure  3.11.  Surviving  Paths  at  tio 
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Figure  3.12.  Surviving  Path  at  time  t 
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Figure  3.13.  Surviving  Path:  Message  Decoded 
3.4  Error  Performance 

As  mentioned  in  Section  IV,  the  atmospheric  channel  is  not  memoryless.  In 
fact,  “the  atmospheric  optical  channel  hcis,  however,  a  very  long  memory  and  a  fade 
can  cause  an  abnormally  large  number  of  erors  that  span  thousands  of  consecutive 
received  channel  bits  (12:993).”  Interleaving  is  a  common  technique  for  overcoming 
the  bursty  errors  that  are  caused  by  channel  memory.  Interleaving  is  a  technique 
where  the  data  bits  are  not  transmitted  successively.  The  first  bit  of  a  byte  is  sent, 

followed  by  the  first  bit  of  the  second  byte,  until  the  first  bit  of  the  byte  is  sent. 

N  is  defined  as  the  interleaving  depth.  After  the  first  bit  of  n  bytes  is  sent,  the 
second  bit  is  sent,  until  all  bits  are  sent.  This  is  illustrated  in  Figure  3.14.  After 
interleaving  to  a  sufficient  depth,  a  burst  of  errors  will  no  longer  effect  more  than 
one  bit  during  the  coherence  time  of  the  channel,  thus  making  the  channel  appear 
memoryless.  Unfortunately,  as  Davidson  points  out: 


Since  the  duration  of  the  fades  is  random,  no  single  maximum  inter¬ 
leaving  depth  can  be  used  to  render  the  channel  completely  memoryless. 
Furthermore,  at  high  source  data  rates,  interleaving  depths  that  corre¬ 
spond  to  time  separations  of  1  ms  between  successive  bits  of  a  codeword 
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Figure  3.14.  Interleaving  a  Data  Stream 

require  the  encoder  and  the  decoder  to  contain  very  large  amounts  of 
memory  (12:993). 

Obviously  the  tradeoff  is  interleaving  depth  for  data  rate.  The  specific  interleaving 
depth  and  data  rate  cissumed  in  this  thesis  are  discussed  in  Chapter  IV. 

As  with  any  digital  communications  link,  when  a  convolutional  coding  scheme 
is  used,  the  probability  of  bit  error  (Pb)  needs  to  be  specified.  Pb  is  a  function  of 
the  code  used,  the  decoding  method,  and  the  probability  of  error  in  the  uncoded 
channel.  The  tool  used  to  account  for  all  these  factors  is  the  transfer,  or  generating 
function  for  the  convolutional  code,  T(D,N).  T(D,N)  is  obtained  from  a  modified 
version  of  the  state  diagram  given  in  Figure  3.5.  This  modified  state  diagram  is 
shown  in  Figure  3.15.  The  a  or  00  state,  can  be  considered  a  beginning  state  and  an 
ending  state,  is  split  up  into  two  states,  a  and  e,  to  explicitly  show  this.  The  state 
transition  paths  are  labeled  with  the  output  codewords  and  a  function  of  D,  where 
the  power  of  D  is  the  Hamming  distance  between  that  codeword  and  the  all  zero 
codeword.  According  to  Sklar, 
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Figure  3.15.  Modified  State  Diagram 


Because  a  convolutional  code  is  a  group  or  linear  code  (28),  there  is  no 
loss  in  generality  in  simply  finding  the  minimum  distance  between  each 
of  the  codeword  sequences  and  the  all-zeros  sequence  (40:338). 

N  is  a  placeholder  that  keeps  track  of  the  number  of  transitions  that  were  caused  by 
an  input  bit  of  1.  The  power  of  N  is  increased  by  one  each  time  a  state  transition 
is  caused  by  a  1.  The  state  equations  for  this  modified  state  diagram  can  now  be 
written  as: 


^6  =  ND'^X^^NXc 
Xc  =  DXi,-VDXi 
Xd  =  NDXi,-^NDXd 

Xe  =  D^Xc  (3.1) 


These  state  equations  can  be  solved  fairly  easily  in  this  case  since  the  number  of 
state  equations  is  given  by  which  is  4  in  this  Ccise  where  K  =  3.  If  >  4  it  is 
more  effective  to  use  matrix  methods  to  solve  the  state  equations.  Two  approaches 


have  been  taken,  differing  only  in  terminology,  not  in  logic.  Sklar  chose  to  leave  Xa 
in  the  equations  and  solve  for  XcjXa  as  the  transfer  function  (40:340).  Clark  and 
Cain  treated  the  input  as  a  unit  impulse,  setting  Xc  =  1  and  then  solved  for  Xc- 
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Using  Clark  and  Cain’s  approach,  the  equations  can  be  written  in  matrix  form  as: 


X  =  AX +Xr 


(3.2) 


where  X  and  Xo  are  2^-1 

X  1  vectors  and  .A  is  a  x  2^  ^  matrix.  For  the  case 
of  the  K  =  Z  encoder  in  Figure  3.2: 

0  //  0  0 
DODO 

A  = 

ND  0  ND  0 
0  1)2  0  0 


and 

’  ND^ 

0 

Xo  = 

0 
0 

Solving  for  X,  X  =  {I  —  A)“’Xo.  The  transfer  function  T(D,N)  is  then  given  by 
X^  in  this  case,  or  by  X^k-i  in  general. 

Viterbi  has  shown  that  Ph  can  be  bounded  using  the  transfer  function  as  shown 
in  equation  3.3  (49:763). 


where  Dq  is  a  function  of  the  uncoded  channel  error  probability.  It  would  be  nice  to 
be  able  to  determine  the  statistics  of  Pb  exactly,  unfortunately  this  can  not  be  done 
since  the  statistics  of  the  path  metrics  can  not  be  evaluated  in  closed  form  (12:996). 
The  following  chapter  will  show  how  these  equations  for  the  bounds  can  be  used  to 
calculate  Pb  curves  for  the  different  coding  schemes. 


3-16 


IV.  Methodology 


Davidson  showed  that  for  a  medium  level  of  turbulence,  5  msec  of  interleaving 
delay  made  the  channel  appear  sufficiently  memoryless.  In  this  case,  the  soft-decision 
encoder’s  measured  bit  errors  closely  agree  with  the  theoretical  curves  for  a  mem¬ 
oryless  channel  (12).  Assuming  a  data  rate  of  10  kbps,  the  required  interleaving 
depth  would  be  50,  whereas  if  a  slower  rate,  like  2  kbps,  could  be  tolerated,  then  the 
required  interleaving  depth  could  be  reduced  to  10.  These  interleaving  depths  can 
be  compared  to  those  that  would  be  required  for  perfect  interleaving.  If  the  time 
scale  of  the  fading  is  cissumed  to  be  1  second,  which  is  typical  for  the  atmosphere, 
then  at  10kbps  the  interleaving  depth  would  have  to  be  1000  and  at  2  kbps  it  would 
have  to  be  200.  These  interleaving  depths  are  impractical  due  to  speed  and  memory 
constraints,  making  Davidson’s  results  extremely  valuable  for  being  able  to  assume 
near  memoryless  conditions  at  substantially  reduced  interleaving  depths. 

Assuming  sufficient  interleaving,  the  generating  function  for  the  encoder  can 
now  be  used  with  the  appropriate  Do  for  predicting  the  bi<^  error  performance.  Using 
the  Poisson  statistics  for  the  photon  arrival  rates,  Davidson  showed  that  for  a  soft- 
decision,  binary  pulse-position  modulation  (BPPM)  channel.  Do  is  computed  to  be 

Do{x)  =  (4.1) 


where  A,  and  Aq  are  the  average  rates  of  arrival  of  the  signal  and  background  photons, 
respectively,  Tq  is  simply  the  bit  duration  time,  and  x  is  the  log-amplitude  of  the 
signal.  In  order  to  compare  the  plots  generated  in  this  thesis  to  Davidson’s  plots, 
this  equation  needs  to  be  written  in  terms  of  Ns  and  No  rather  than  Aj,  Aq  and  To. 
Bringing  2b/2  inside  the  radical,  equation  4.1  can  be  rewritten  as  equation  4.2. 


Do{x)  =  e 


Nse^X+Ng 


(4.2) 
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As  discussed  in  Chapter  II,  scintillation  effects  cause  x  to  be  normally  distributed. 
To  account  for  the  scintillation  effects,  Z(x)  must  be  integrated  over  this  normal 
density  to  obtain  the  final  expression  for  D  given  in  equation  4.3. 


The  integration  in  equation  4.3  was  first  done  numerically  using  MathStation  (32). 
Unfortunately,  the  resulting  Pj’s  were  three  to  four  orders  of  magnitude  higher  than 
those  reported  by  Davidson.  Since  x  is  normally  distributed,  the  first  check  was  to 
ensure  that  the  area  above  and  below  the  numerical  limits  of  integration  were  in  fact 
negligible.  As  the  range  of  integration  was  increased,  the  value  of  Do  decreased,  ver¬ 
ifying  that  the  integration  wais  being  done  incorrectly.  The  version  of  MathStation 
being  used  was  determined  to  be  designed  for  an  older  version  of  the  SUN  operating 
system.  This  seemed  to  a  likely  explanation  for  the  incorrect  integration.  Ko\,ever,  a 
new  version  of  MathStation  was  received  which  was  supposed  to  be  compatible  with 
present  version  of  the  operating  system.  Unfortunately,  this  new  version  gave  the 
same  results  as  the  old  version.  The  final  option,  which  fortunately  was  successful, 
WcLS  to  use  the  same  twenty  point  gauss-hermite  integration  that  Davidson  used.  A 
C  program  was  written  to  perform  the  summation.  It  is  included  in  Appendix  ^ 
This  program  Wcis  checked  on  a  zero  mean,  unit  variance,  normal  distribution.  The 
result  was  within  .001%  of  the  correct  result  of  one.  After  multiplying  the  normal 
distribution  by  Do(x)  ^  given  by  equation  4.2  and  integrating,  a  arithmetic  under¬ 
flow  was  encountered.  Since  the  problem  seemed  to  be  ocurring  in  the  exponential, 
equation  4.3  was  rewritten  using  equation  4.1  with  To  =  1,  instead  of  equation  4.2. 
The  probability  of  bit  error  curves  are  then  obtained  by  substituting  equation  4.3 
into  the  generating  function  and  plotting  over  a  range  of  values  for  Ng  and  spe¬ 
cific  values  of  Nq  and  cr^.The  resulting  Pb's  were  within  one  order  of  magnitude  of 
Davidson’s  results. 
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Chapter  III  showed  the  derivation  of  the  state  equations  for  the  case  of  K  =  3. 
Writing  out  a  state  diagram  for  cases  where  K  >A  gets  extremely  cumbersome.  The 
method  used  in  this  thesis  was  to  write  out  all  possible  states  twice  and  then  put 
either  a  1  or  a  0  in  front,  representing  an  input  bit  of  1  or  0.  The  output  codeword 
was  then  determined  by  the  taps  and  the  contents  of  the  register.  The  method  is 
illustrated  in  equation  4.4. 

000  Xa  Xa 

100  Xe  — ^  Xc 

001  Xb  X, 

101  Xb  — ^  Xc 
010  Xc  Xb 

110  Xc  — ^  Xj 

oil  Xj  Xb 

111  Xj  — >•  Xj 

The  double  arrow  represents  a  transition  caused  by  a  zero.  Using  these  equations, 
the  state  equations  can  be  written  as: 


00 

11 

11 

00 

01 

10 

01 

10 


(4.4) 


Xb  =  DXc  +  BXj 
Xc  =  ND^Xa-f-NXb 
Xd  =  NDXc  +  NBXj 

Xc  =  B^Xb  (4.5) 


These  state  equations  are  the  same  as  those  given  in  Chapter  III  except  for 
the  6  and  c  rows  being  swapped.  This  is  simply  a  result  of  the  b  and  c  states  being 
swapped  from  10  to  01  and  vice  versa. 

The  codes  examined  in  this  thesis  are  derivatives  of  a  code  developed  by  Massey 
and  Costello.  They  developed  an  algorithm  to  determine  generating  polynomials 
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that  would  produce  non-catastrophic,  nonsystematic,  optimal  codes.  One  of  the 
codes  they  discovered  has  a  constraint  length  of  48  and  is  represented  octally  by 
533,533,676,737,355,3  and  733,533,676,737,355,3.  Fortunately,  they  also  showed  that 
shorter  codes  with  these  same  desirable  characteristics  could  be  obtained  by  simply 
dropping  elements  off  of  the  right  end  until  the  desired  constraint  length  is  reached 
(30).  Interestingly,  and  conveniently,  the  K  =  3  code  obtained  in  this  manner  is 
the  same  as  the  example  mentioned  in  Chapter  III.  This  provided  an  excellent 
means  for  verifying  the  results  obtained,  as  they  could  be  checked  against  the  work 
done  by  Davidson.  The  method  could  then  be  extended  to  the  longer  constraint 
codes  with  reasonable  confidence  as  to  the  accuracy  of  the  results.  The  extension 
is  straightforward  conceptually,  however,  the  number  of  calculations  required  for  an 
exact  solution  becomes  immense.  In  order  to  obtain  an  exact  solution,  the  matrix 
(I-A)  must  be  inverted  symbolically  in  terms  of  D  and  N.  This  was  done  on  a  Sun  3 
workstation  using  Mathematica  (31),  a  software  package  that  can  handle  symbolic 
mathematics.  For  the  K  <5  case,  the  inversion  was  done  with  no  appreciable  time 
delay.  Unfortunately,  for  K  =  6,  after  a  week,  the  computer  had  not  finished  the 
inversion.  Obviously,  a  numerical,  approximate  solution  was  needed.  Viterbi  states 
in  his  article  on  convolutional  codes  that 


...  a  generating  function  calculation  is  basically  just  a  matrix  inversion, 
which  can  be  performed  numerically  for  a  given  value  of  D.  The  derivative 
at  N  =  1  can  be  upper  bounded  by  evaluating  the  first  difference 

lT{D,l  +  c)-T{D,l)]/e  (4.6) 


for  small  e  (49:761) 


The  Pfc  curve  in  Figure  4.1  was  produced  using  Pb  <  for  K  =  3, 

where,  —  2Do)^.  Figure  4.2  shows  the  curve  produced 

by  the  the  numerical  method  with  e  =  .01.  These  two  curves  are  overlaid  in  Fig¬ 
ure  4.3.  The  actual  data  that  was  used  to  generate  the  plots  is  shown  in  Ta- 
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Figure  4.1.  Pi,  Curves  from  Exact  T(D,N) 


J.  !« !M! !!  M  irpiJI  !!!* ! !  rrnti  |t  |  In  it  itf  ni  1 1  li  1 1  ji  li !!  nlH !!  I  HI !!!!  iTl! !!!?!!!!  111!  H  !M 

0.001  l!!!p(i!!!!n!!!!!)!!!!!i!!!!!U!!!!!i!!!!!ji!!!!!(!!!!l!j!!!!!!j!!!!!j!i!!U!!!!!U!!!!n 

0 , 0001  fiMfi  ! ! rS !l !■!•*!** !t " !  !>•! it !!«!!! i.  11  !l !! II !* 'I i>> M >!•{}'"■!«! II IS!! !!*!•} MSfl 

. \mZlZlZ^  Z  m  mm  m  m^ZlZZZZ^m  •Z  •l. 

^  *  . . \ . . . . . . . --  I . y . , . . . . . •**.». 

le-08  II !!!! jil!!! !!!!!!!! nHiniUniu^!!!!! !!!!!!!!!!!!!!! Hi! !!!!!!!!!! If!! !!!!}!!!!!! 

le-09  r-7"i-"V""'i:''"i . . f-'-r'— 

1  2  3  4  5  6  7  8  9  10  11  12  13  14  15 

Ns 


Figure  4.2.  Pb  Curves  from  Approximate  Method 
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Table  4.1.  Comparison  of  Exact  and  Numerical  Methods  for  Finding 


Ns 

%  Error 

1 

1.219574370237228  IQ-^ 

1.257684126095466  10"^ 

3.10 

2 

5.408651034564673  lO""* 

5.440526112330407  lO"-* 

0.58 

3 

2.068339667884289  10"® 

2.073866818209250  10"® 

0.26 

4 

1.625352476271095  10"® 

1.627815938605986  10"® 

0.15 

5 

2.010669397505719  10“^ 

2.012617133662639  lO”'^ 

0.09 

6 

3.339719190100789  10"® 

3.341940128714635  10"® 

0.06 

7 

6.676329368437486  10"® 

6.679512188932890  10"® 

0.05 

8 

1.507209044081476  10"® 

1.507738692119392  10“® 

0.04 

9 

3.725211292170222  10”^° 

3.726196035074627  lO"^® 

0.03 

9.956548639013677  lO”^^ 

9.958562668240980  10"“ 

0.02 

11 

2.867396086188684  10“^^ 

2.867847057830152  10"” 

0.02 

12 

8.891537688029467  lO''^^ 

8.892641949738258  lO-'"^ 

0.01 

13 

2.965991554500035  lO'^^ 

2.966286851056265  lO-^^ 

0.01 

14 

1.061641638028914  lO-^^ 

1.061727606677203  lO-^^ 

0.01 

15 

1.648835809525375  10“^® 

1.648927667019529  10"^® 

0.01 

ble  4.1.  The  percent  error  using  the  numerical  approximation,  was  calculated  as 
^exact:^2isen«2ijl00%.  It  Can  be  seen  in  the  plots  and  from  the  data  that  the  numer¬ 
ical  approximation  is  a  good  one  for  Ns  >  4.  For  iV*  <  4  neither  method  provides 
realistic  results  as  P^’s  greater  than  one  and  even  less  than  zero  were  obtained. 
Davidson’s  plots  show  the  curves  asymptotically  approaching  1  for  Ns  ~  2.  He 
makes  no  mention  of  the  fact  that  the  numbers  do  not  make  practical  sense  for  this 
region  of  low  SNR. 

This  method  can  now  be  used  to  calculate  the  P^  curves  for  all  constraint 
lengths  <  7.  The  D  parameter  from  equation  4.3  Wcis  calculated  for  Ns  ranging 
from  1  to  15.  The  D  parameter  was  then  used  in  a  MATRIXx  (33)  program  that 
evaluated  equation  4.6.  The  MATRIXx  command  files  that  were  written  are  shown 
in  Appendix  A.  For  K  —  8,  the  huge  number  of  calculations  and  resulting  roundoff 
errors  made  the  matrix  appear  singular  within  MATRIXx’s  precision.  However, 
based  on  the  results  obtained  from  K  =  6  and  K  =  7,  as  shown  in  Figures  5.1 
through  5.9,  the  curves  start  getting  closer  and  closer  with  increeising  K.  Since  7  is  a 
common  constraint  length,  it  seemed  like  a  good  point  to  stop  based  on  diminishing 
returns. 


4-7 


V.  Results 


The  results  of  Chapter  IV’s  calculations  are  presented  here.  Tables  of  the  raw 
data  are  given  in  the  appendices,  while  the  Pb  curves,  are  presented  in  this  chapter. 

The  MATRIXx  command  files  from  Appendix  A  were  used  to  analyze  codes 
with  constraint  lengths  from  3  to  7.  The  D  parameter  that  was  used  as  an  input 
was  calculated  for  no  turbulence,  {<jfnx  =  .01),  low  turbulence,  =  .2)  and  for 
medium  turbulence,  —  -6),  as  well  cis  for  three  different  background  light  levels. 
No  =  .1,  1,  and  10  photons/bit.  parameter  Davidson  uses  to  describe  the 

strength  of  the  turbulence.  He  defines  X{t)  as  the  received  signal,  by: 

(5.1) 

He  defines  A*  as  the  “information  bearing  intensity  modulation  impressed  on  the 
laser  beam  at  the  transmitter  (12:994).”  af^x  then  defined  cis: 

=  (5-2) 

Using  equation  5.1  and  cissuming  that  Xs{t)  is  deterministic,  equation  5.2  can  be 
rewritten  as: 


alx  =  Var{lnX{t)} 

=  Var  {lnXs{t)  +  2x{t)} 

= 

=  4a^  (5.3) 

Tables  of  the  D  parameter  are  given  in  Appendix  C.  Tables  of  the  resulting  Pb  data 
are  given  in  Appendix  D.  Although  the  calculations  were  made  in  each  case  for  Ng 
ranging  from  1  to  15,  in  each  case  there  was  a  certain  N^^,^  below  which  the  signal  to 


5-1 


noise  ratio  Wcis  too  low  for  this  method  to  work.  Below  Pb  could  go  negative 

or  even  above  1.  Data  points  for  Ns  <  Ns^,„  were  not  plotted.  They  were  retained 
in  the  tables  of  data  found  in  Appendix  D.  Ns„,„  was  determined  for  each  case  as 
the  point  where  the  curve  broke  continuity  with  decreasing  Ns- 

When  the  results  of  these  calculations  are  comp'  '•d  to  Davidson’s  results,  there 
is  a  difference  of  w  one  half  of  an  order  of  magnitude,  i  he  cause  for  this  difference  is 
the  evaluation  of  equation  4.3.  The  evaluation  of  equation  4.3  is  very  dependent  on 
the  method  of  integration,  and  the  precision  of  the  calculator.  The  goal  of  this  thesis 
was  to  provide  the  communications  engineer  a  means  of  determining  the  coding  gain 
available  from  longer  constraint  length  codes.  The  plots  in  this  chapter  provide  that. 
The  absolute  value  of  Pt  is  not  known  anyway  since  these  are  only  bounds  on  the 
error.  The  atmosphere  is  very  difficult  to  model  with  any  degree  of  reliability.  The 
ability  of  these  curves  to  accurately  represent  the  correct  value  of  Pb  is  dependent 
upon  the  level  of  turbulence,  which  is  me^lsured  with  cr^.  is  directly  proportional 
to  C^,  which  can  vary  3  to  4  orders  of  magnitude. 

Another  factor  that  is  somewhat  hidden  in  this  analysis,  is  the  dependence  of 
Pb  on  path  length.  This  dependence  is  accounted  for  in  cr^.  Throughout  this  thesis 
a  certain  <t^  has  been  assumed,  cr^  can  be  calculated  from  C^,  (9:36)  cis: 

<tI  =  (5.4) 

where  k  is  still  the  wavenumber,  which  can  be  given  cis  a>/c  or  2ir/A,  and  L  is  the 
pathlength.  Unfortunately,  equa'ion  5.4’s  validity  breaks  down  as  the  pathlength 
increases  or  C^,  increeises.  Davidson  alluded  to  this  problem  in  his  discussion  of  the 
data  obtained  for  .6  (12:999). 

Even  with  these  limitations  a  lot  of  information  is  shown  in  these  plots.  Fig¬ 
ures  5.?  through  5.9  show  that  Pb  can  be  varied  over  5  orders  of  magnitude  with  the 
use  of  5  constraint  lengths  and  a  fairly  small  range  of  Ns-  For  example.  Figure  5.4 
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shows  that  for  Ns  ranging  from  5  to  10,  a  dynamic  range  of  only  3dB,  Pb  can  be 
varied  from  5  •  10“^  with  K  =  3  at  iV,  =  5  to  less  than  10“^  with  K  =  7  and  Ng  =  10. 
Figures  5.3,  5.6  and  5.9  show  that  if  the  background  light  gets  too  high,  turbulence 
does  not  matter  anymore.  Inordinately  large  values  of  Ng  are  required  to  get  Pb 
low  enough  to  be  of  any  use.  This  is  shown  for  constraint  lengths  of  3,5,  and  7  in 
Figures  5.11  and  5.12  respectively. 

Another  important  consideration  is  demonstrated  in  these  plots.  When  select¬ 
ing  codes,  the  last  shift  register  should  be  tapped.  Otherwise,  the  encoder  does  not 
produce  a  maximal  length  code  and  the  Pb  is  the  same  as  that  obtained  from  a  code 
of  the  same  length  tis  the  length  from  the  first  register  in  the  code  to  the  last  register 
that  h2is  a  tap.  This  is  shown  in  Figures  5.1  through  5.9.  The  Pb  curves  for  K  =  Z 
is  the  same  as  K  =  A  and  the  curves  for  K  =  6  overlays  K  =  7.  Seven  is  a  common 
constraint  length  for  a  convolutional  code,  unfortunately  for  this  particular  code, 
there  is  no  performance  gain  over  a  constraint  length  of  six. 
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K  =3  - 

K  =4 - 

K  =5  ---- 

K  =6  . 

K  =7 - 

Figure  5.1.  P6  for  K  =  3  thru  7  With  No  Turbulence,  =  .1 
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K  =  3  - 

K  =  4 - 

K  =  5  ---- 

K  =  6  . 

K  =  7  -- 

Figure  5.2.  Pi,  for  K  =  3  thru  7  With  No  Turbulence,  No  =  1 
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Figure  5.5.  Pb  for  K  =  3  thru  7  With  Low  Turbulence,  No  =  l 
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K  =  3  - 

K  =  3 - 

K  =  3  ---- 

K  =  3  . 

K  =  3 - 

Figure  5.6.  P6  for  K  =  3  thru  7  With  Low  Turbulence,  No  =  10 
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K  =  7 - 

Figure  5.8.  Pj  for  K  =  3  thru  7  With  Medium  Turbulence,  No  =  1 
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K  =  6  . 

K  =  7  — • 

Figure  5.9.  Pb  for  K  =  3  thru  7  With  Medium  Turbulence,  No  =  10 
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Figure  5.10.  Comparing  Turbulence  and  Background  Effects  for  K  =  3 
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Figure  5.11.  Comparing  Turbulence  and  Background  Effects  for  K  =  5 
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Figure  5.12.  Comparing  Turbulence  and  Background  Effects  for  K  =  7 
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VI.  Conclusions 


The  objective  of  this  thesis  was  to  show  that  adaptive  ED  AC  could  combat  the 
time  varying  fading  caused  by  atmospheric  turbulence.  The  Pi,  plots  in  Chapter  V 
show  that  this  is  definitely  acheivable.  With  a  small  number  of  codes  with  constraint 
lengths  varying  from  3  to  7  and  with  a  laser  transmitter  with  an  intensity  range  of 
«  6dB  (3  <  <  12),  Pi,  can  be  varied  over  a  range  of  5-6  orders  of  magnitude.  The 

communications  engineer  can  now  assume  a  value  of  and  a  path  length,  calculate 
<T^  and  pick  the  constraint  length,  /V,  combination  that  will  give  the  required  Pi,. 

This  thesis  was  the  first  attempt  at  analyzing  the  clear-air  turbulent  atmo¬ 
spheric  optical  channel  for  adaptive  ED  AC  applications.  There  are  numerous  oppor¬ 
tunities  for  further  research.  An  experimental  setup  similar  to  Davidson’s  would  be 
invaluable  for  verifying  the  predictions  made  in  this  thesis.  It  would  enable  analysis 
of  the  effects  of  varying  and  path  length.  Better  algorithms  for  analyzing  longer 
constraint  lengths  could  be  developed  and  other  codes  could  then  be  investigated. 
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Appendix  A.  MATRIXx  Programs 


The  following  sections  contain  the  programs  that  were  used  in  MATRIXx  to 
calculate  the  Pt,  for  K  —  Z  through  K  =  7.  The  main  program  is  the  same  for  all 
cases  except  for  calling  functions  that  depend  on  the  constraint  length.  The  main 
program  for  K  =  3  is  shown  in  section  A.l.  The  only  change  to  the  main  program 
for  a  different  constraint  length  is  changing  the  number  in  the  T101K3  and  T100K3 
functions. 

A.l  Main  Program 

d  =[.356998  .245867  .165898  .111127  .0742772  .0496584  .0332469  .0223035  ... 
.0149948  .0101029  .00682063  .00461301  ,00312477  .00211944  .00143909]; 
long  e; 

for  i  =  1 : 15, . . . 

dpar  =  d(i) .... 

n  =  1.01, .. . 

define  'tl01k3>,. . . 

tlOl  =  tl01k3(n,dpar) , . . . 

n  =  1.0, . . . 

define  'tl00k3’,,.. 

tlOO  =  tl00k3(n,dpar) , . . . 

pb(i)  =  (tlOl  -  tl00)/.01, 

end; 

print  ('pbk3s0’ ,PB); 
return 

A. 2  K3  Functions 

//TlOl  =  T101k3(N,DPAR) 

A  =  Cl,-N,0,0;-DPAR,1,-DPAR,0;-N+DPAR,0,1-N*DPAR,0;0,-DPAR**2,0,1] ; 

X  =  [N*DPAR**2,0,0,03; 

T  =  INV(A)*X’; 
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TlOl  =  T(4) ; 

RETF 

//TlOO  =  T100k3(N,DPAR) 

A  =  [1 , -N , 0 , 0 : -DPAR , 1 , -DPAR, 0 ; -K*DPAR , 0 , 1-N*DPAR, 0 ; 0 , -DPAR**2 , 0 , l3 
X  =  [M*DPAR**2, 0,0,0]; 

T  =  mV(A)*X'; 

TlOO  =  T(4); 

RETF 

A. 8  K4  Functions 

//TlOl  =  T101k4(N,DPAR) 

A  =  [1,-DPAR**2,-DPAR**2, 0,0, 0,0,0 
0,1,0, -DPAR, -DPAR,0, 0,0 
0,0, 1,0,0, -DPAR, -DPAR, 0 
-N*DPAR**2, 0,0, 1,0, 0,0,0 
0,-N,-N, 0,1, 0,0,0 
0,0,0, -N*DPAR , -N*DPAR ,1,0,0 
0 , 0 , 0 , 0 , 0 , -N*DPAR , 1-K*DPAR , 0 
-1,0,0,0,0,0,0,13; 

X  =  [0,0,0,N*DPAR**2, 0,0, 0,03; 

T  =  INV(A)*X'; 

TlOl  =  T(8); 

RETF 

//TlOO  =  T100k4(N,DPAR) 

A  =  Cl,-DPAR**2,-DPAR**2,0,0,0,0,0 
0,1,0, -DPAR, -DPAR,0, 0,0 
0,0, 1,0,0, -DPAR, -DPAR, 0 
-JI*DPAR’f*2, 0,0, 1,0, 0,0,0 
0,-N,-N, 0,1, 0,0,0 
0,0,0,-N*DPAR,-K*DPAR, 1,0,0 
0 , 0 , 0 , 0 , 0 , -K*DPAR, 1-N*DPAR , 0 
-1,0, 0,0, 0,0, 0,1]; 
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X  =  [0,0,0,N*DPAR**2,0,0,0,0]; 

T  =  INV(A)*X»; 

T100=  T(8): 

RETF 

A. 4  K5  Functions 

//TlOl  =  T.01k5(N,DPAR) 

A  =  [1,-1,-DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,1,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0 
0,0,1,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0 
0,0,0,1,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0,0 
0,0,0,0,1,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0,0 
0,0,0,0,0,1,0,0,0,0,0,-DPAR,-DPAR, 0,0,0 
0,0,0,0,0,0,1,0,0,0,0,0,0,-DPAR,-DPAR,0 
-N, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0 
0,-N*DPAR**2,-N, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0,0 
0,0,0,-N,-N*DPAR**2, 0,0, 0,0, 1,0, 0,0, 0,0,0 
0,0,0,0,0,-N,-N*DPAR**2, 0,0, 0,1, 0,0, 0,0,0 
0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,0, 1,0, 0,0,0 
0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,1, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 1,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,l-N*DPAR,0 
-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,1]; 

X  =  [0,0,0,0,0,0,0,M*DPAR**2,0,0,0,0,0,0,0,03; 
T  =  INV(A)*X'; 

TlOl  =  T(16); 

RETF 


//TlOO  =  T100k5(K,DPAR) 

A  =  Cl,-1,-DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,1,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0 
0,0,1,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0 
0,0,0,1,0,0,0,-DPAR,-DPAR,0,0,0,0,0,0,0 
0,0,0,0,1,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0,0 


A-3 


0,0,0,0,0,1,0,0,0,0,0,-DPAR,-DPAR, 0,0,0 
0,0,0,0,0,0,1,0,0,0,0,0,0,-DPAR,-DPAR,0 
-N, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0 
0,-N*DPAR**2,-N, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0,0 
0,0,0,-N,-N*DPAR+*2, 0,0, 0,0, 1,0, 0,0, 0,0,0 
0,0,0,0,0,-H,-N*DPAR»*2, 0,0, 0,1, 0,0, 0,0,0 
0,0,0,0,0,0,0,-K*DPAR,-N*DPAR, 0,0, 1,0, 0,0,0 
0,0,0,0,0,0,0,0.0,-K>i‘DPAR,-H*DPAR, 0,1, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 1,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,1-N+DPAR,0 
-N+DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,1]; 

X  =  CO,0,0,0,0,0,0,N*DPAR**2,0,0,0,0,0,0,0,0]; 
T  =  INV(A)*X’: 

TlOO  =  T(16); 

RETF 


A. 5  K6  Functions 


//TlOl  =  T101k6(N,DPAR) 

A  =  Cl,-DPAR**2, -1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0 
0,1,0,-1,-DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,1,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
0,0,0,1,0,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
0,0,0,0,1 ,0,0, 0, 0,-1, -DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,1,0,0,0,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,(  ,0,0, 0,0 
0,0,0,0,0,0,1,0,0,0,0,0,0,-1,-DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,0,0, -DP AR,-DPAR, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR,0 
-N, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 


A-4 


0,-N,-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,-N»DPAR**2,-H, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,-K,-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,-N,-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,-K*DPAR**2,-N, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,-N,-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR**2,-N, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-H*DPAR,-N*DPAR, 0,0, 0,0, 1,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR,0,0,  ,1,0, 0,0, 0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,0, 1,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,C,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,1, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 1,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,l-N*DPAR,0 
-DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,1] ; 

X  =  C0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,N*DPAR**2, 0,0, 0,0, 0,0, 0,0,0, 0,0, 0,0, 0,0,0]; 
T  =  INV(A)*X'; 

TlOl  =  T(32); 


//TlOl  =  T101k6(N,DPAR) 

A  =  [1,-DPAR**2, -1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0 
0,l,0,-l,-DPAR+*2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,1,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
0,0,0,1,0,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
0,0,0,0,l,0,0,0,0,-l,-DPAR+*2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,1,0,0,0,0,0,-DPAR**2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
0,0,0,0,0,0,l,0,0,0,0,0,0,-l,-DPAR*+2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,-DPAR,-D?AR, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0, 0,0,0 


A-5 


0,0,0,0,0,0,0,0,0,0,0,0,1.0,0,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR, 0,0, 0,0,0 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,0, 0,0, 0,0, 0,0, -DP AR,-DPAR, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-DPAR,-DPAR,0 
-N, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,-N,-N*DPAR**2,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0  '',0,0, 0,0, 0,0, 0,0, 0,0, 0,0 
0,0,0,-K*DPAR**2,-N,0,0,0,0,0,0,0,0,0,0,0,0,1,0,U,0,0,0,0,0,0,0,0,0,0,0,0 
0,0,0,0,0,-N,-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,-N,-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,-N*DPAR**2,-N, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,-N,-N*DPAR**2, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR*’t'2,-N,  0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-M*DPAR, 0,0, 0,0, 1,0, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,0, 0,1, 0,0, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAP  ‘',0,1, 0,0, 0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 0,1, 0,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,-N*DPAR, 1,0,0 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-N*DPAR,1-N*DPAR,0 
-DP AR**2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,1] ; 

X  =  CO,0,0,0,0,0,0,0,0,0,0,0,0,0,0,N*DPAR**2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] ; 
T  =  INV(A)+X’; 

TlOl  =  T(32): 

RETF 


A. 6  K7  Functions 

//TlOl  =  T101k7(N,DPAR) 

Z  =  0*0NES(8,8); 

A  =  [1,-DPAR**2,-DPAR**2,0,0,<  ,0,0 
0,1  0.-DPAR**2,-DPAR**2, 0,0,0 
'.,0.  .0,-1, -1,0 
0,0,C,_,j,0,0,-l 


0,0, 0,0, 1,0, 0,0 
0, 0.0, 0,0, 1,0,0 


A-6 


0,0, 0,0, 0,0, 1,0 
0,0,0,0,0,0,0,13; 

B  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-1,0, 0,0, 0,0, 0,0 
0 , -DPAR**2 , -DPAR**2 , 0 , 0 , 0 , 0 , 0 
0, 0,0, -DP AR**2,-DPAR**2, 0,0,0 
0,0, 0,0, 0,-1, -1,0 
0,0,0,0,0,0,0,-DPAR**2] ; 

C  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR**2, 0,0, 0,0, 0,0,0]  ; 

D  =  [0,-1, -1,0, 0,0, 0,0 
0,0, 0,-1, -1,0, 0,0 
0 , 0 , 0 , 0 , 0 , -DPAR**2 , -DPAR**2 , 0 
0,0,0,0,0,0,0,-DPAR**2 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 ; 

E  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR**2, 0,0, 0,0, 0,0,0 
0,-1, -1,0, 0,0, 0,0 
0,0, 0,-1, -1,0, 0,0 
0,0,0,0,0,-DPAR**2,-DPAR**2,0 
0,0,0,0,0,0,0,-DPAP.3; 


A-7 


0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR, 0,0, 0,0, 0,0,0] ; 

G  =  Co, -DPAR, -DPAR, 0,0, 0,0,0 
0,0,0, -DPAR, -DPAR,0, 0,0 
0, 0,0, 0,0, -DPAR, -DPAR, 0 
0,0, 0,0, 0,0,0, -DPAR 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03: 

H  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR, 0,0, 0,0, 0,0,0 
0, -DPAR, -DPAR, 0,0, 0,0,0 
0,0,0, -DPAR, -DPAR, 0,0,0 
0, 0,0, 0,0, -DPAR, -DPAR, 0 
0,0, 0,0, 0,0,0, -DPAR]; 

I  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0.0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR, 0,0, 0,0, 0,0,0] : 

J  =  [0,0, 0,0, 0,0, 0,0 

0,0, 0,0, 0,0, 0,0 


A-8 


0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N*DPAR**2, 0,0, 0,0, 0,0,0] ; 

K  =  [0,-DPAR,-DPAR,0,0,0,0,0 
0,0,0,-DPAR,-DPAR, 0,0,0 
0,0,0,0,0,-DPAR,-DPAR,0 
0,0,0,0,0,0,0,-DPAR 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 : 

L  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR, 0,0, 0,0, 0,0,0 
0,-DPAR,-DPAR, 0,0, 0,0,0 
0,0,0, -DPAR, -DPAR, 0,0,0 
0, 0,0, 0,0, -DPAR, -DPAR, 0 
0,0,0,0,0,0,0,03 ; 

M  =  [0,-N,-N, 0,0, 0,0,0 
0,0,0,-N,~N, 0,0,0 
0 , 0 , 0 , 0 , 0 , -NfDPAR**2 , -N*DPAR**2 , 0 
0,0,0,0,0,0,0,-N*DPAR**2 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 ; 

NM  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N*DPAR**2, 0,0, 0,0, 0,0,0 


A-9 


0,-N,-N, 0,0, 0,0,0 
0,0,0,-N,-N, 0,0,0 
0 , 0 , 0 , 0 , 0 , -N*DPAR**2 , -N*DPAR**2 , 0 
0,0,0,0,0,0,0,-N3; 

0  =  Co, 0,0, 0,0, 0,0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N, 0,0, 0,0, 0,0,0]; 

P  =  [0,-N*DPAR**2,-N*DPAR**2, 0,0, 0,0,0 
0,0,0, -N*DPAR**2 , -N*DPAR**2 ,0,0,0 
0,0,0,0,0,-N,-N,0 
0,0,0,0,0,0,0,-N 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 : 

Q  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N, 0,0, 0,0, 0,0,0 

0 , -N*DPAR**2 , -N*DPAR**2 , 0 , 0 , 0 , 0 ,0 
0,0,0, -N*DPAR**2 ,-N*DPAR**2 ,0,0,0 
0,0,0,0,0,-N,-N,0 
0,0,0,0,0,0,0,-N*DPAR3; 

R  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
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0,0, 0,0, 0,0, 0,0 
-N*DPAR,0,0,0,0,0,0,03; 

S  =  C0,-N*DPAR,-N*DPAR,0,0,0,0,0 
0,0,0,-N*DPAR,-N*DPAR, 0,0,0 
0 , 0 , 0 , 0 , 0 , -N+DPAR, -N*DPAR , 0 
0,0,0,0,0,0,0,-N*DPAR 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 
T  =  Co, 0,0, 0,0, 0,0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N*DPAR, 0,0, 0,0, 0,0,0 
0 , -N*DP AR , -N+DPAR , 0 , 0 , 0 , 0 , 0 
0,0,0, -N+DPAR , -N*DPAR ,0,0,0 
0,0,0,0,0,-N*DPAR,-N*DPAR,0 
0,0,0,0,0,0,0,-N*DPAR3 ; 

U  =  Cl, 0,0, 0,0, 0,0,0 
0,1, 0,0, 0,0, 0,0 
0,0, 1,0, 0,0, 0,0 
0,0,0, 1,0, 0,0,0 
0,0, 0,0, 1,0, 0,0 
0,0, 0,0, 0,1, 0,0 
0,0, 0,0, 0,0, 1,0 
-N*DPAR,0,0,0,0,0,0,13; 

V  =  C0,-N*DPAR,-N*DPAR,0,0,0,0,0 
0,0,0,-N*DPAR,-N*DPAR, 0,0,0 
0 , 0 , 0 , 0 , 0 , -NfDPAR , -H*DP AR , 0 
0,0,0,0,0,0,0,-N*DPAR 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 : 
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W  =  Cl, 0,0, 0.0. 0,0,0 
0,1, 0,0, 0,0, 0,0 
0,0, 1,0, 0,0, 0,0 
-N*DPAR,0,0,1,0,0,0,0 
0 , -N*DP AR , -N+DPAR,  0 , 1 , 0 , 0 , 0 
0,0,0, -N*DPAR . -N+DPAR, 1,0,0 
0 , 0 , 0 , 0 , 0 , -N*DP AR , 1-N*DP AR , 0 
0,0, 0,0, 0,0, 0,1] : 

GEN  =  [A,B,C,Z,Z,Z,Z,Z 

Z,EYE(8),D,E.F,Z,Z,Z 

Z,Z,EYE(8),Z,G,H,I.Z 

J,Z,Z,EYE(8),Z,Z,K,L 

H,NH,0,Z,EYE(8),Z,Z,Z 

Z,Z,P,Q,R,EYE(8),Z,Z 

Z,Z,Z,Z,S,T,U,Z 

0,Z,Z,Z,Z,Z,V.W]; 

XI  =  [0.0,0,0,0.0,0,03; 

X2  =  [0,0,0.0.0.0,0,03: 

X3  =  [0,0,0,0,0,0,0,03: 

X4  =  [0,0,0,0,0,0,0,N*DPAR**2]: 
X5  =  [0,0,0,0,0,0,0,03: 

X6  =  [0,0,0,0,0,0,0,03; 

X7  =  [0,0,0,0,0,0,0,03; 

X8  =  [0,0,0,0,0,0,0,03; 

X  =  [X1,X2,X3,X4,X5,X6,X7,X8]: 
T  =  INV(GEN)*X'; 

TlOl  =  T(64); 
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//TlOO  =  T100k7(N.DPAR) 

Z  =  0*0NES(8,8); 

A  =  [1,-DPAR**2,-DPAR**2, 0,0, 0,0,0 
0,1,0, -DPAR**2 , -DPAR**2 ,0,0,0 
0,0, 1,0, 0,-1, -1,0 
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0,0,0,1|0,0,0,-1 
0,0, 0,0, 1,0, 0,0 
0,0, 0,0, 0,1, 0,0 
0,0, 0,0, 0,0, 1,0 
0,0, 0,0, 0,0, 0,1] ; 

B  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-1,0, 0,0, 0,0, 0,0 
0 , -DPAR**2 , -DPAR**2 , 0 , 0 , 0 , 0 , 0 
0,0,0, -DPAR**2 , -DPAR**2 ,0,0,0 
0,0, 0,0, 0,-1, -1,0 
0,0,0,0,0,0,0,-DPAR**2]: 

C  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR**2, 0,0, 0,0, 0,0,0] ; 

D  =  [0,-1, -1,0, 0,0, 0,0 
0,0, 0,-1, -1,0, 0,0 
0 , 0 , 0 , 0 , 0 , -DPAR**2 , -DPAR**2 , 0 
0,0,0,0,0,0,0,-DPAR**2 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0] ; 

E  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0. 0,0, 0,0 
0,0, 0,0 

-DPAR-  0,0, 0,0,0, 0,0 
0,-1, -1,0, 0,0, 0,0 
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0,0, 0,-1, -1,0, 0,0 
0 , 0 , 0 , 0 , 0 , -DPAR**2 , -DPAR**2 , 0 
0,0,0,0,0,0,0,-DPAR3 ; 

F  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR, 0,0, 0,0, 0,0,03; 

G  =  [0, -DPAR, -DPAR, 0,0, 0,0,0 
0,0,0, -DPAR, -DPAR, 0,0,0 
0, 0,0, 0,0, -DPAR  ..-DPAR, 0 
0,0, 0,0, 0,0,0, -DPAR 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 : 

H  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR, 0,0, 0,0, 0,0,0 
0,-DPAR,-DPAR, 0,0, 0,0,0 
0,0,0, -DPAR, -DPAR,0, 0,0 
0,0, 0,0,0, -DPAR, -DPAR,0 
0,0,0,0,0,0,0,-DPAR3 ; 

I  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0. 0.0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
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-DPAR, 0,0, 0,0, 0,0,0] ; 

J  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N*DPARt>*2, 0,0, 0,0, 0,0,0]  ; 

K  =  [0, -DPAR, -DPAR, 0,0, 0,0,0 
0,0,0, -DPAR, -DPAR,0, 0,0 
0, 0,0, 0,0, -DPAR, -DPAR, 0 
0,0, 0,0, 0,0,0, -DPAR 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0] ; 

L  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-DPAR, 0,0, 0,0, 0,0,0 
0, -DPAR, -DPAR, 0,0, 0,0,0 
0,0,0, -DPAR, -DPAR,0, 0,0 
0, 0,0, 0,0, -DPAR, -DPAR, 0 
0,0, 0,0, 0,0, 0,0] : 

M  =  [0,-N,-N, 0,0, 0,0,0 
0,0,0,-N,-N, 0,0,0 
0 , 0 , 0 , 0 , 0 , -N*DPAR**2 , -N*DPAR**2 , 0 
0,0,0,0,0,0,0,-N*DPAR**2 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0] ; 

NM  =  [0,0, 0,0, 0,0, 0,0 
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0,0, 0,0, 0,0, 0,0 

0,0, 0,0, 0,0, 0,0 

-N*DPAR**2, 0,0, 0,0, 0,0,0 

0,-N,-N, 0,0, 0,0,0 

0,0,0,-N,-N, 0,0,0 

0 , 0 , 0 , 0 , 0 ,  -Ht>DPAR**2 ,  -H*DPAR**2 .0 

0,0,0,0,0,0,0.-N]; 

0  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N, 0,0, 0,0, 0,0,0] ; 

P  =  [0,-K*DPAR**2,-N*DPAR**2, 0,0, 0,0,0 
0,0,0,  -N*DPAR*4‘2 ,  -N*DPAR*t2 ,0,0,0 
0,0,0,0,0,-N,-K,0 
0,0,0,0,0,0,0,-N 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0] ; 

Q  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N, 0,0, 0,0, 0,0,0 

0 , -N*DPAR**2 , -N*DPAR**2 , 0 , 0 , 0 , 0 , 0 
0,0,0, -N*DPAR**2 , -N*DPAR**2 .0,0,0 
0,0,0,0,0,-H,-N,0 
0,0,0,0,0,0,0,-N*DPAR] ; 

R  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
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0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N*DPAR, 0,0, 0,0, 0,0,0] ; 

S  =  C0,-N*DPAR,-H*DPAR,0,0,0,0,0 
0,0,0, -K*DPAR , -H*DPAR ,0,0,0 
0 , 0 , 0 , 0 , 0 , -N*DP AR , -N*DP AR , 0 
0,0,0,0,0,0,0,-N*DPAR 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0,0,0,0,0,03 ; 

T  =  [0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
-N*DPAR, 0,0, 0,0, 0,0,0 
0,-NfDPAR,-N*DPAR, 0,0, 0,0,0 
0,0,0, -NfDPAR , -N*DPAR ,0,0,0 
0,0,0,0,0,-N*DPAR,-N*DPAR,0 
0,0,0,0,0,0,0,-N*DPAR]; 

U  =  Cl, 0,0, 0,0, 0,0,0 
0,1, 0,0, 0,0, 0,0 
0,0, 1,0, 0,0, 0,0 
0,0, 0,1, 0,0, 0,0 
0,0, 0,0, 1,0, 0,0 
0,0, 0,0, 0,1, 0,0 
0,0, 0,0, 0,0, 1,0 
-N*DPAR,0,0,0,0,0,0,1]; 

V  =  [0,-N*DPAR,-N*DPAR,0,0,0,0,0 
0,0,0,-N*DPAR,-N*DPAR, 0,0,0 
0 , 0 , 0 , 0 , 0 ,-N*DPAR, -N*DPAR, 0 
■0,0,0,0,0,0,0,-II*DPAR 
0,0, 0,0, 0,0, 0,0 
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0,0, 0,0, 0,0, 0,0 
0,0, 0,0, 0,0, 0,0 
0,0,0, 0,0, 0,0,0]; 

W  =  Cl, 0,0, 0,0, 0,0,0 
0,1, 0,0, 0,0, 0,0 
0,0, 1.0, 0,0, 0,0 
-H*DPAR, 0,0, 1,0, 0,0,0 
0 , -N*DPAR , -NKDP AR , 0 , 1 , 0 , 0 , 0 
0 , 0 , 0 , -N*DP AR , -N*DP AR , 1 , 0 , 0 
0 , 0 , 0 , 0 , 0 , -N*DPAR, 1-N+DPAR ,0 
0,0,0,0,0,0,0,11 : 

GEN  =  LA,B,C,Z,Z,Z,Z,Z 

Z,EYE(8),D,E,F,Z,Z,Z 

Z,Z,EYE(8),Z,G,H,I,Z 

J,Z,Z,EYE(8),Z,Z,K,L 

H,NM,0,Z,EyE(8),Z,Z,Z 

Z,Z,P,Q,R,EYE(8),Z,Z 

Z,Z,Z,Z,S,T,U,Z 

o,z,z,z,z,z,v,w]: 

XI  =  [0.0,0,0,0,0,0,01; 

X2  =  [0,0,0,0,0,0,0.01; 

X3  =  [0,0,0,0,0,0,0,01; 

X4  -  [0,0,0,0,0,0,O.N*DPAR**2l; 
X5  =  [0,0,0,0,0,0.0,01; 

X6  =  iC. 0,0, 0,0, 0,0,01; 

X7  =  [0,0  0,0,0,0,0.01; 

X8  =  [0,0,0,0,0,0,0,01; 

X  =  [X1,X2,X3,X4,X5,X6,X7,X8]; 
T  --  IMV(GEN)*X’; 

TlOO  =  T(64); 
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Appendix  B.  C  Programs 


#include  <stdio.h> 
#include  <math.h> 


^define  pi  3.1415926535898 
Mafine  t  1.0 


double  xvecC]  =  { 

-5,3874808900112» 

-4.6036824495507, 

-3.9447640401156, 

-3.3478545673832, 

-2.7888060584281, 

-2.2549740020893, 

-1.7385377121166, 

-1.2340762153953, 

-0.737^.  *^7285454, 

-0.2453407083009, 

0.2453407083009, 

0.7374737285454, 

1.23407621539^ i, 

1.7385377121166, 

2.2549740020893, 

2.7888060584281, 

3.3478545673832, 

3.9447640401156, 

4.6036824495507, 

5.3874808900112 

>: 


double  wvecC3  =  { 
0.8985919614532, 


0.7043329611769 


0.6222786961914, 

0.5752624428525, 

0.5448517423644, 

0.5240803509486, 

0.5096790271175, 

0.4999208713363, 

0.4938433852721, 

0.4903215006667, 

0.4909215006667, 

0.4938433852721, 

0.4999208713363, 

0.5096790271175, 

0.5240803509486, 

0.5448517423644, 

0.5752624428525, 

0.6222786961914, 

0.7043329611769, 

0.8985919614532 

}; 


double  myint(ns,  no,  sx) 
int  ns; 

double  no,  sx; 

•C 

int  i; 

double  Is  =ns/t; 
double  g,  sum  =  0.0; 

for  (i  =  0;  i  <  20;  i  ++) 

{ 

g  -  exp(-((pow((sqrt((ls  +  l)*exp(2*xvec[i3 )  +  no) 

-  sqrt(no)),2.0))/(t/2,0)))*exp  (-(xvecCi]  *  xvecCi3)/(2  *  sx)) 
/  sqrt  (2  *  pi  *sx); 
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sum  +=  wvecCi]  ♦  g; 

> 

return  stun; 

} 


void  make_matrixx_table(no,  sx,  k, filename) 
double  no,  sx; 
int  k; 

char  ♦filename; 

i 

int  ns; 

FILE  *fout; 

font  =  fopen  (filename, 
fprintf  (fout,  **d=C"); 
for  (ns  =  0;ns<14;ns++) 

fprintf  (fout,  •'y,16.14e  ..An",  myint  (ns,  no,  sx)); 
fprintf  (fout,  "y,16.14e]  ;\n",  myint  (15,  no,  sx)); 
fprintf  (fout,  "long  e;\nfor  i  =  1: 15, . . .\ndpar 
=  d(i) , . . .\nn  =  1.01, . . .\n") ; 
fprintf  (fout,  "define  'tl01ky,d' , . .  .\ntl01  = 
tl01ky,d(n,dpar) , . .  .\n",  k,  k); 

fprintf  (fout,  "n  =  1.0, . .  .\ndefine  'tl00ky.d\ . .  .\n" ,  k); 
fprintf  (fout,  "tlOO  =  tl00ky.d(n,dpar) , . .  .\npb(i)  = 

(tlOl  -  tl00)/.01,\n",  k); 

fprintf  (fout,  "end;\nprint  (^pby,s\pb) ;\n",  filename); 


f close  (fout); 

return; 

> 
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void  main  () 

i 


make^matrixx^table  (0.1, 
make.matrixx^table  (1.0, 
make^matrixx.table  (10.0, 
make^matrixx^table  (0.1, 
make^matrixx^table  (1.0, 
make^matrixx^table  (10.0, 
make^matrixx.table  (0.1, 
make^matr ixx_table  (1.0, 
make_matr ixx_t able  (10.0, 
make.matrixx^table  (0.1, 
make_matr ixx.table  (1.0, 
ma)ce_matrixx.table  (10.0, 
make^matrixx^table  (0.1, 
make^matrixx^table  (1.0, 
make^matrixx^table  (10.0, 
make^matrixx^table  (0.1. 
make^matrixx^table  vl.O, 
make^matrixx^table  (10.0, 
make^matrixx^table  (0.1, 
make^matrixx^table  (1.0, 
make^matrixx.table  (10.0, 
make^matrixx^table  (0.1, 
make.matrixx^table  (1.0, 
make^matr ixx^t able  (10.0, 
make^matrixx^table  (0.1, 
make^matr ixx^t able  (1.0, 
make_matrixx_table  (10.0 
make^matrixx^table  (0.1, 
make^matrixx^table  (1.0, 
make^matrixx^table  (10.0 


.01,  3,  ’'kSsO"); 

.01,  3,  "kSsONol"); 
.01,  3,  "k3sONolO"); 
.2,  3,  '‘k3s2*'); 

.2,  3,  "k3s2Nol»): 

.2,  3,  "k3s2NolO"); 
.6,  3,  "kSse’O; 

.6,  3,  ’•k3s6Nol*'); 

.6,  3,  ’*k3s6Nol0'‘); 
.01,  4,  ••k4s0"); 

.01,  4,  ’‘k4sONol“); 
.01,  4,  "k4sONolO"); 
.2,  4,  »k4s2‘'); 

.2,  4,  "k4s2Nol"); 

.2,  4,  '•k4s2NolO’'); 
.6,  4,  "k4s6'‘); 

.6,  4,  ••k4s6Nol'*); 

.6,  4,  "k4s6NolO’'); 
.01,  6,  "kSsO"); 

.01,  6,  "kSsONol*'); 

.  .01,  5,  "kSsONolO'O; 
.2,  5,  "k5s2‘*); 

.2,  5,  "k5s2Nol"); 

,  .2,  5,  "k5s2NolO"); 
.6,  5,  "k5s6**); 

.6,  5,  "k5s6Nol"); 

,  .6,  5,  “kSseNolO") ; 
.01,  6,  "k6s0"); 

.01,  6,  "kesONoi'O; 

,  .01,  6,  "kesONolO"); 
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make^matrixx^table  (0.1,  .2,  6,  *'k6s2"); 
make^matrixx^table  (1.0,  .2,  6,  **k6s2Nol'*) ; 
make^matrixx^table  (10.0,  .2,  6,  "k6s2NolO"); 
make^matrixx^table  (0.1,  .6,  6,  **k6s6"); 
make^matrixx^table  (1.0,  .6,  6,  *'k6s6Nol")  5 
make_matrixx_table  (10.0,  .6,  6,  "k6s6NolO**) ; 
nake.matrixx^table  (0.1,  .01,  7,  "kTsO") ; 
inake_inatrixx_oable  (1.0,  .01,  7,  *'k7s0Nol"); 
wake^matrixx^table  (10.0,  .01,  7,  "k7sONolO") 


make.matrixx^table  (0.1, 

.2, 

7. 

"k7s2"); 

make.inatrixx.tablo  (1.0, 

.2. 

7. 

"k7s2Nol"); 

make_matrixx_table  (10.0, 

.2 

.  7, 

''k7s2NolO": 

make.matrixx.table  (0.1, 

.6, 

7. 

"k7s6''): 

make^matrixx^table  (1.0, 

.6, 

7. 

"k7s6Nol"): 

make^matrixx^table  (10.0,  .6,  7,  *'k7s6NolO'') ; 
exit  (0); 

} 
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Appendix  C.  D  Parameter  Data 


Tab 


e  C.l.  D  Parameter  Data 


No  Tiirb. 

Low  Tlirb. 

Medium  Tbrb. 

No 

No 

No 

.1 

1 

10 

.1 

1 

10 

.1 

1 

10 

6.838  c-2 

1.279  e-l 

1.801  e-1 

3.759  e-1 

6.232  e-l 

8.773  e-1 

4.080  e-1 

5.794  e-1 

7.818  e-1 

2.493  e.2 

6.954  c-2 

1.509  e-1 

1.847  e-1 

3.851  e-1 

7.181  e-1 

2.655  e-1 

4.180  e-1 

6.441  e-1 

9.618  e-3 

3.948  e-2 

1.199  e-1 

1.052  e-1 

2.585  e-1 

5.861  e-1 

1.950  e-1 

3.293  e-1 

5.531  e-1 

3.643  e-3 

2.244  e-2 

9.272  e-2 

6.572  e-2 

1.814  e-1 

4.804  e-1 

1.526  c-1 

2.704  e-1 

4.836  e-1 

1.346  e-3 

1.231  e-2 

7.280  e-2 

4.411  e-2 

1.309  e-1 

4.003  e-1 

1.243  e-1 

2.281  e-1 

4.304  e-1 

4.871  e-4 

6.496  e-3 

5.793  c-2 

3.115  e-2 

9.757  e-2 

3.395  e-1 

1.040  e-1 

1.967  e-1 

3.891  e-1 

1.734  c-4 

3.305  e-3 

4.636  e-2 

2.274  e-2 

7.511  e-2 

2.912  e-1 

8.882  e-2 

1.726  e-1 

3.555  e-1 

8 

6.097  c-5 

1.630  e-3 

3.694  e-2 

1.696  e-2 

5.948  e-2 

2.510  e-1 

7.692  e-2 

1.534  e-1 

3.268  e-1 

9 

2.120  c-5 

7.837  e-4 

2.910  e-2 

1.287  e-2 

4.810  e-2 

2.166  e-1 

6.744  e-2 

1.377  e-1 

3.015  e-1 

10 

7.309  e-6 

3.682  e-4 

2.258  e-2 

9.911  e-3 

3.944  e-2 

1.872  e-1 

5.976  e-2 

1.245  e-1 

2.792  e-1 

11 

2.499  e.6 

1 .696  e-4 

1.723  e-2 

7.740  e-3 

3.263  e-2 

1.622  e-1 

5.344  c-2 

1.131  e-1 

2.596  e-1 

12 

8.494  «-7 

7.679  e-5 

1.293  e-2 

6.132  e-3 

2.714  e-2 

1.412  e-1 

4.816  e-2 

1.034  e-1 

2.423  e-1 

13 

2.869  e-7 

3.422  e-5 

9.554  e-3 

4.928  c-3 

2.267  c-2 

1.236  e-1 

4.370  e-2 

9.490  c-2 

2.272  e-1 

14 

9.641  e.8 

1.504  e-5  1 

6.948  e-3 

4.0)6  e-3 

1.900  e-2 

1.091  e-1 

3.987  e-2 

8.750  c-2 

2.139  e-1 

15 

1.074  e.8 

2.800  e-6 

3.518  e-3 

2.769  e-3 

1.352  e-2 

8.702  e-2 

3.365  e-2 

7.531  c-2 

1.918  c-1 
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Appendix  D.  data 


The  following  five  tables  document  the  Pb  results.  These  tables  include  all 
calculated  values,  even  if  Ng  <  .  Each  table  contains  nine  data  files.  Each  table 

contains  data  files  for  one  constraint  length,  three  different  levels  of  turbulence,  and 
three  different  levels  of  background  light.  The  data  from  these  files  is  what  was  used 
to  generate  the  plots  in  Chapter  V. 
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Table  D.l.  Pb  Data  for  K 


Pb  For  if  =  3 


a 

1.76118240665752610-“ 

4.904566565602125  10-“ 

3.396456792650220 10-» 

2 

1.017354924451575  10-® 

1.921734348660993  10-® 

1.224958709911291 10-* 

3 

8.397319078211453 10-** 

1.048272598181580  10-® 

3.331418923555024 10-® 

4 

6.472600105900685 10-*® 

5.969356340095553  10"® 

8.678203923716110 10“® 

5 

4.433856655295366 10-*® 

2.910710890111120 10-*® 

2.438359965463750 10“® 

6 

2.747571891027707 10-*® 

1.173029370202268  10-** 

7.469991627962507 10-® 

7 

1.571882969699214  10-*® 

3.971810008318049  lO'*® 

2.379987415479274  IQ-® 

8 

8.431511677911093 10-®® 

1.157133576249314  10“** 

7.469081871613898 10-® 

9 

4.288591615981567 10-®* 

2.962219528746474  10“*® 

2.223642918456624  10-® 

10 

2.085955427953614  10-®® 

6.778258264562654  lO"*® 

6.159539365369669 10-® 

11 

9.764961200283160  10“®® 

1.405706140023401  10-*® 

1.575958127837387 10-® 

12 

4.421784567919114  10-®* 

2.672098384729986  10“®* 

3.720772305832319 10-’® 

13 

1.944613366786925 10“®® 

4.699436067252185  10“®® 

8.121053675096788 10-** 

14 

8  332955515699094  10“®® 

7.706607805107916  10-®® 

1.643361047373747  10"** 

15 

1.43030787733953210-*® 

1.723468515013948  10-®« 

5.430471929799078  lO"'® 

8.120055014739710 10-' 
3.931892282434991 10“^ 
1.70052967012164710-^ 
1.433177452290168 10-® 
1.84564816827370510-^ 
3.142249121934885  10-® 
6.384467711497623 10-® 
1.457612075765014  10-* 
3.631626773978695  10-*« 

10  9.763100487750492  10-^* 

11  2.823742398542454  10-** 

12  8.784030800346874  10-*^ 

13  2.937118961913594 10“** 

14  1.053208729634933 10-*2 

15  1.639785267277624  10-*® 


2.379629443127217  10 
3.992748811609001  10-® 
5.441960393269649 10-* 
1.306331513774349 10-* 
4.196162484604749  10-® 
1.606196044084013 10-® 
6.916003649504374  10“® 
8  1  3.250888820356039 10-® 
1.639654355303673 10-® 

10  8.77117151662408710-’’ 

11  4.931641502445706 10-^ 

12  2.892425784560671  10-^ 

13  1.75816834875907510-’’ 

14  1.10157433991161710—^ 

15  4.654814508057136 10“® 


8.487659687208049  10-' 
1.082286983237514  10- * 
3.314796936897471  10-® 
3.537424615151867  10-* 
5.576841310259403  10-® 
1.137747267107984  10-® 
2.870491703555186  10-® 
8.557950975777767  10-^ 
2.873272713027236  lO”^ 
1.042681984068186  10"^ 
3.972916049777352  10“® 
1.562163872009560  10-® 
6.285513800467057  10“® 
2.581357298837393  10-® 
4.655390295906908  10-*® 


1.195674781507047  10 
3.546305080024785  10-* 
2.154722825671706  10-® 
4.556196067477102  10-® 
1.441293220643628  10"® 
5.726198895502450  10-^ 
2.642157249132219  10-* 
1.347803248022064  10-* 
7.362405206986415  10-® 
4.224100177599823  10-® 
2.517739630895961  10-® 
1.549766482689044  10-® 
9.818900461017002  10-® 
6.3S95S065256972G  10—® 
2.911160472784485  10-® 


1.677902410456666  10“ 
7.696636284990427 10-* 
1.10680133136736010® 
6.342054048982953 10* 
1.801994950094393 10-* 
2.848829344308801 10-® 
7.875677543651548  10-® 
2.709052136780571  10-® 
1.038846812090549  10-® 
4.265466560153942  10-* 
1.850424482966563  IQ-* 
8.457867363303756  10-® 
4.074988505470772  10-® 
2.069959431916303 10“® 
6.209172645633017  10“® 


9.271490168948904  10 
7.867306842671201  10"* 
1.795235665863260 10® 
2.259715112525575  10® 
5.964349212153219 10-* 
1.234850134312394 10"* 
4.454524327495591  10“® 
2.015718277251066 10"® 
1.032264965290606 10“® 
5.745679712175315 10“® 
3.412336555122953 10“® 
2.139584527041704  10“® 
1.405331785242145  10-® 
9.603416451709336 10—- 
4.925816172545805 10-< 
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Table  D.2.  Ph  Data  for  K 


Pb  For  A  =  3 


1.761182406657526  lO-" 
1.017354924451575  10"® 
8.397319078211453 10”" 
6.472600105900685  10“* *® 
4.433856655295366 10“*® 
2.747571891027707 10“*^ 
1.571882969699214  10"*- 
8.43151167791109310-*® 
4.288591615981567  lO"*^ 
2.085955427953614 10"*® 
9.764961200283160  lO"*® 
4.421784567919114  10”®* 
1.944613366786925  10"®® 
8.332955515699094  10“®® 
1.430307877339532  lO'^® 


4.904566565602125  10““ 
1.921734348660993 10“® 
1.048272598181580  lO"* 
5.969356340095553  10“® 
2.910710890111120  10“'® 
1.173029370202268 10-” 
3.971810008318049  10“‘® 
1.157133576249314  10-*< 
2.962219528746474  10“*® 
6.778258264562654  10"*® 
1.405706140023401  10“*® 
2.672098384729986  lO”** 
4.699436067252185  10”*® 
7.70660780S10791G  lO"*® 
1.723468515013948  10“*® 


3.396456792650220 10“’ 
1.224958709911291  10“* 
3.331418923555024  10“® 
8.678203923716110 10“® 
2.438359965463750  10“® 
7.469991627962507 10“* 
2.379987415479274  10“* 
7.469081871613898 10“® 
2.223642918456624  10“® 
6.159539365369669 10“® 
1.575958127837387 10“® 
3.720772305832319 10“*® 
8.121053675096788 10““ 
1.643361047373747  10““ 
5.430471929799078  10“*® 


8.120055014739710  lO"'* 
3.931892282434991  10-< 
1.700529670121647  10“® 
1.433177452290168  ;r.-® 
1.845643168273705  10“* 
3.142249121934885  10“® 
6.384467711497623  10"® 
1.457612075765014  10“® 
3.631626773978695  10“*® 
9.763100487750492  10““ 
2.823742398542454  10“” 
8.784030300346874  10"** 
2.937118961913594  10"** 
1.053208729634933 10“** 
1.639785267277624  10“*® 


5.487659687208049  10“' 
1.082286983237514  lO"' 
3.314796936897471  10“® 
3.537424615151867  10“^ 
5.576841310259403  10“® 
1.137747267107984  10“® 
2.870491703555186  10“® 
8.557950975777767  10“* 
2.873272713027236  10”* 
1.042681984068186 10”* 
3.972916049777352  10“® 
1.562163872009560  10“® 
6.285513800467057  10”® 
2.581357298837393  10“® 
4.655390295906908 10“*® 


1.677902410456666 10” 
7.696636284990427  10“* 
1.106801331367360 10® 

•  6.34205404898295310* 
1.801994950094393  10“* 

2.848829344308801  10“* 
7.875677543651548  10“® 
2.709052136780571 10“® 
1.038846812090549  10“® 
4.265466560153942 10“^ 
1.850424482966563 10-< 
8.457867363303756  10“® 
4.074988505470772 10“® 
2.069959431916303 10“® 
6.209172645633017 10“® 


1  I  2.379629443127217 1 
3.992748811609001 1 
5.441960393269649 1 
1.306331513774349  1 
4.196162484604749 1 
1.606196044084013 1 
6.916003649504374 1 
3.250888820356039 1 
9  1.639654356303673 1 

10  8.7711715166240871 

11  4.9316415024457061 

12  2.892425784560671 ] 

13  1.758168348759075 1 

14  1.101574339911617] 

15  4.6548145080571361 


1.195674781507047  10” 
3.546305080024785  10“* 
2.154722825671706  10“* 
4.556196067477102  10"® 
1.441293220643628 10“® 
5.726198895502450 10-< 
2.642157249132219  lO"* 
1.347803248022064  10“^ 
7.362405206986415  10“® 
4.224100177599823  10“® 
2.517739630895961  10“® 
1.549766482689044  10“® 
9.818900461017002  10“® 
6.389580652569720 10“® 
2.911160472784485  10“® 


9.271490168948904  10“' 
7.867306842671201  10“* 
1.79523566586326010® 
2.259715112525575 10* 
5.964349212153219 10“* 
1.234850134312394 10“* 
4.454524327495591  10“* 
2.015718277251066  10“* 
1.032264965290606  10“* 
5.745679712175315  10“® 
3.412336555122953  10“® 
2.139584527041704  10“® 
1.405331785242145 10“® 
9.603416451709336  10-‘ 
4.925816172545805  lO'^ 


Table  D.3.  Pb  Data  for  K  =  5 


Pb  For  /C  =  3 


BS 

No  IWb.,  No  =  *1 

No  Torb.,  No  =  1 

No  Tlirb.,  No  =  10 

3.5888451767737?.';  10^^ 

1.301184503768433  10“’ 

1.680465130267246 10“' 

7.6690387613870S:  10-‘® 

3.980163380726879  10“’ 

5.207015185079078 10“’ 

2.4464194364364J5  lO-*^ 

1.247034326894963  10“® 

1.147238600864936 10“’ 

7.146243354404256 10“** 

4.051058721145679  10“*® 

2.36651262V595794  10“® 

1.808640834230535 10~'^ 

1.085829575287451  10“** 

5.278314017989475  10“’ 

4.056102188233781  lO"*® 

2.308813460411382  10“*’ 

1.295504847913651 10“’ 

8.262857754515191 10“®® 

3  977720907063443  10“*’ 

3.317402180471925 10“* 

8 

1.557876481989643 10“®® 

5.717712803349629  10“*’ 

8.318987758782499 1C“® 

9 

2.755690890.547344  10“®® 

7.034889943650598  10“*® 

1.954607725392788  10“® 

10 

4.619803437705312  10“’* 

7.563678831881658  10“®* 

4.206064372755250 10“*® 

11 

7.400280816438396 10“’* 

7.226038558594962  10“®’ 

8.219595965743579 10“** 

12 

1.137425173407888 10“’« 

6.218221367756666  10“®’ 

1.457409248334728 10“** 

13 

1.716024625870967  10“” 

4.874214330825844  10“®’ 

2.350232400630777  10“*® 

14 

1.131821909637099 10“*® 

3.513042449737500  10“®® 

3.459435483874404  10“*’ 

15 

1.336613602433585  10“*® 

1.462674764641170  10”^^ 

5.789118843527627  10“*’ 

Ns 

Low  Turb.,  No  =  1 

Low  T\irb.,  No  =  10 

n 

1.11250574743286710“* 

3.645008647716608  10“® 

-1.297001939504244  10“* 

2 

1.986239552879131 10“* 

1.735913277627548  10“* 

-1.34348651995505610"® 

3 

5.220252950621978 10“® 

2.218733031621265  10“’ 

6.048894486972989 10“® 

4 

2.810136834317699 10“’ 

1.760390946552966  10“* 

4.23216404254579710“* 

5 

2.449585741767166 10“® 

2.091636336584438  10“’ 

4.313034274341834  10“* 

6 

2.955676813040052  10“® 

3.255309287777114  10“® 

2.719520709096677 10“® 

7 

4.390605358812497 10“*® 

6.403854195633279  10“’ 

5.929377005524828 10“’ 

8 

7.485702453305598 10“** 

1.522872728180790  10“’ 

1.766789109960466  10“’ 

9 

1.415655429742433 10“** 

4.152646824989154  10“® 

5.988339013774150  10“* 

10 

2.930919049102203 10“*® 

1.239097453498596  10“® 

2.179569747221685  10“* 

11 

6.621426663135924 10“*’ 

3.912452551063778  10“® 

8.376239766728474  10“’ 

12 

1.632017401396818 10“*’ 

1.281217910795444  10“® 

3.392294739359754  10“’ 

13 

4.385698529831749 10“** 

4.309259054064914  10“*® 

1.451365962390224  10“’ 

14 

1.281530040667071  10“** 

1.484465853269795  10“*® 

6.573972458374902  10“« 

15 

1.376239930903162  10“*’ 

1.906482908117359  10“** 

1.594551718750755  10“« 

Ns 

7.885245478859597  10“* 

6.61264873459757610“^ 

.5.279573991508191  10“^ 

2.738044105650970 10“’ 

2.446917658423788  10® 

2.518265434996725 10“® 

2.876358773413516 10“* 

1.930808588654748 10“® 

9.278919042043851  10“® 

5.608346337297015 10“’ 

3.179356740378413  10“’ 

3.788577707523860  10“* 

5 

1.501815469693323 10“’ 

8.669529860511201  10“* 

7.888675059683534  10* 

6 

4.881993564268271 10“® 

3.047802679370277  10“* 

2.158718824316941  10“* 

7 

1.810875825513764  10“® 

1.260725060765081  10“* 

4.815927308105591  10“® 

8 

7.420511338777283 10“’ 

5.814292437280116  10“’ 

1.781422226651555  10“® 

9 

3.296925255066699 10“’ 

2.888256612107963  10“’ 

8.097102819865998  10“’ 

10 

1 .568091169499945 10“’ 

1.513479966289268  10“’ 

4.136945599605393 10“’ 

11 

7.903780365572642  10“® 

8.268861750717344  10“« 

2.292611438951678  10“’ 

12 

4.185664232452487 10“® 

4.681223928538208 10“® 

1.353884242526594  10“’ 

13 

2.31171211357377310“® 

2.737054316171640 10“® 

8.425694484550083 10“* 

14 

1.323110279446483  10“® 

1.649374424006569  10“® 

6.478898389952230 10“  • 

15 

4.728023384668811  10“® 

6.511444169874516  10“’ 

2.56799453975670.3 10“* 
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Table  D.4.  Pb  Data  for  K  = 


Pb  For  K  =  3 


2.382127307176335  lO"® 
1.885098299156822 10“” 
2.338866274725782  10-*< 
2.597894551151169 10“‘^ 
2.432694545137668 10-^® 
1.975465431803489  lO"*® 
1.432550075619665 10-^® 
9.427133212499651  10"®® 
4.390171200035802  10"®® 
-6.829907668649090 10“®® 
-6.102907131704052  10“®® 
-3.190332457977668 10-<* 
-1.600907898305692  10-‘< 
-7.746733446803838 10-<® 
-1.650446056129215  10"®* 
Low  "Ilirb.,  No  =  .1 
5.743664664500233  lO"^ 
3.606910679826663 10“® 
5.309820419209465  lO”’^ 
1.793677628647456  10“® 
1.056452915164334  10“® 
9.050546551583628  10““ 
9.852770990751123  10“*® 
1.257210509321701  lO”*® 
1.807980926719250 10“*® 
2.886941955968257 10“*^ 
5.100324452423429 10“*® 
9.969320132346569  10“*® 
2.154608109469527 10"*® 
5.133228135621588  lO"*® 
3.805204741211403  lO"*® 


8.241289142822711  lO"* 
7.538455212590924  10“^ 
5.540986245500433 10“® 
8.324699564395067 10"® 
1.807291092696138 10“® 
4.913101530064253 10“® 
1.556237210469711  10"® 
5.531498954225660 10“® 
2.158870038080265 10“® 
9.116049391667563 10“® 
4.116355274622969  10-® 
1.968021342102399 10“® 
9.877141033860166  lO"*® 
5.165878504579067 10“*® 
1.562538235892088 10“*® 


2.231174284297292  10-“ 
2.685947826201847  10“® 
4.822446476683164  10"*® 
8.071983476719887 10"*® 
1.327454414624348 10“*® 
1.493787718307394  10"*® 
1.311940624735994 10"*® 
9.314232260044080 10"®® 
5.510861355375329 10"®® 
2.785015473498183 10“®* 
1.226271236584166 10"®® 
4.773036242179096 10“®® 
1.470028954599524  10"®* 
8.757901655161567 10"®* 
-1.351955786795744  10"®® 


2.659093328550022  lO"-' 
9.973576791457144  10“® 
5.909969593869129  10-‘ 
3.135830529243702  10“® 
2.653360716068365  10“® 
3.071416476355582  10“® 
4.662676338210941  10“® 
8.811736495408450  10“® 
1.950008395977639  10“® 
4.786807601157403 10"*® 
1.253868662464033  lO"*® 
3.424219821258997 10“" 
9.640711118410.  41  lO”*® 
2.789452413443463 10"*® 
2.557150877844377  10"*® 


3.576630460481109  lO"-' 
8.798089953412937  10® 
7.231698424961953  10"® 
8.961280180452228  10-‘ 
1.990548074458836  10-< 
5.925685001413348  10“® 
2.129316359944198  10"® 
8.680138704528980  10"® 
3.857084205125528  10”® 
1.823572134293640  10”® 
9.049358413713859  10"® 
4.679675146589373  10“® 
2.512122523783557  10“® 
1.396505996971246  10“® 
4.753641391946674  lO"® 


2.971355535087142  10"“ 
7.638074861784263  10"® 
1.322674451991810  10“® 
2.122301074413344  10"® 
3.726342876442348  10"® 
7.304223399876416  10“® 
1.502435890064750  10“® 
3.013036943800496  10“*® 
5.595893982260899  10““ 
9.372634491312656  10"*® 
1.401822796484459  10“*® 
1.870238706283164  10“*® 
2.232070504786458  10-*< 
2.393246054009976  10“*® 
2.032193798223647  10"*® 


1.270522175501332  10 
1.498327504435987 10“® 
3.397367088680131  10"® 
1.320314549901232  10“* 
3.332709571871491  10”* 
1.077408259820985  10“® 
1.842621245179735  10"® 
4.540243831457460  10-< 
1.296463118689234  10“^ 
4.017051017898418  10“® 
1.325418499208964  10"® 
4.647741271967976  10“® 
1.737108236745757  10"® 
6.936649046796719  10“® 
1.342335958718455  10“® 


1.046862481187821  lO"-' 
2.342114954357391  10“® 
4.427587837375485  10"® 
1.218374264515466  10"® 
6.013387701980970 10® 
1.316470239336894  10"* 
2.105792179801268  10“® 
6.586058016557492  lO"® 
2.642124493447580 10“® 
1.215162366372452  10“® 
6.138316514995245  10“^ 
3.335956822324857  10"^ 
1.925825601366948  10"^ 
1.169515771512415  10"^ 
4.858330376914143  10"® 


Table  D.5.  Pb  Data  for  K 


Pb  For  K  =  3 


2.382127344268520  lO"® 
1.885098299347114  lO"" 
2.338866263838570  10-*< 
2.597894751999962  10“*^ 
2.432695279028801  10“*® 
1.975484127927358  lO"*® 
1.433290331900348  10"®® 
9.499167544345708 10"®® 
5.843646755618153 10“®® 
3.376597671210300  10"®® 
1.849243461916648  10“®® 
9.667029172946557  10-‘® 
4.850911495913263 10"^® 
2.347337778672049  lO"*® 
5.001016612943'132  10"®® 


5.824923614211278  10-“ 
3.606989252995464  10"® 
5.309822229150763 10“^ 
1.793677649692370  10"® 
1.056452915912844  10”® 
9.050546551417740  10-“ 
9.852770990838970  10““ 
1.257210508153758 10““ 
1.807980925588390  lO"*® 
2.886941951223575  10““ 
5.100324448461079  10““ 
9.969320140927946  10"*® 
2.154608172258122  lO"*® 
5.133227699673478  10-“ 
3.805204902308587  10-*® 


8.900673734152466 10-' 
-  M31603395345 10-< 

6.541169474610012  10“® 
8.324742695118008 10-® 
1.807293166444064  10-® 
4.913103080562563 10-^ 
1.556237367233361  IQ-^ 
5.531499153427659 10-® 
2.158870068548048 10-® 
9.116949446076360  10-® 
4.116355285741972  10-® 
1.968021344677371  10-® 
9.877141040231267  10-J® 
3.165873506375393 10—*® 
1.562538236124325  10->® 


2.231177437699840  10"® 
2.685947873331373 10-® 
4.822446478409413  10-'® 
8.971983475273210  10-*® 
1.327454412475624  10“'® 
1.493787728667977  10-'® 
1.311940683553941  10-"^ 
9.314237940559646 10-®® 
5.510852355491087  IQ-®® 
2.784751489577215  10-®‘ 
1.225746588087130 10-®® 
4.775306103028195  10-®® 
1.668395817368469  10-®' 
5.285132649300303  10-®' 
4.096565130858280  10-®® 


2.417387030075810  10-“ 
1.018044453478983  IQ-' 
5.911832549244163  10-' 
3.135890128443548  10-® 
2.653365166996777 10-® 
3.071417084369847  10-® 
4.662676479826909  10-® 
8.811736546414640  10-® 
1.950008398448117  10-® 
4.786807602464710  10-'® 
1.253868662576982  10-'® 
3.42421982159130710-“ 
9.640711115665297  10-'® 
2.789452413741822  10-'® 
2.557150879581083  10-'® 


3.313992078286046  10-“ 
1.202220882925985  lO' 
7.252539989458179  10-® 
8.965427753358741  10-' 
1.990773337652379  IQ-' 
5.925894130770415  10-® 
2.129344086599585  IQ-® 
8.680185565183570  10-® 
3.857093567532652  10-® 
1.823574245373506 10-® 
9.049363647592486 10-® 
4.679676553764520  10-® 
2.512122931061629 10-® 
1.3365061 232871 58  10—® 
4.753641539174235  10-® 


8.478781911559763  10-® 
2.105949429850110 10-® 
4.137056498324937 10-® 
1.147267622247223 10-' 
4.703434192633561 10® 
1.349380162431932 10-' 
2.120017548430369 10-® 
6.603621975059663 10-® 
2.645351729816206 10-® 
1.215904752701084  10-® 
6.140320737379708 10-' 
3.336572619354721  10-' 
1.926036783228483 10-' 
1.103595335934243 1C-' 
4.858471781701328  10-® 
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